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ABSTRACT 

Satellite galaxies in groups and clusters are more likely to have low star formation 
rates (SFR) and lie on the 'red-sequence' than central ('field') galaxies. Using galaxy 
group/cluster catalogs from the Sloan Digital Sky Survey Data Release 7, together 
with a high-resolution, cosmological TV-body simulation to track satellite orbits, we 
examine the star formation histories and quenching timescales of satellites of M sta r > 
5 x 10 9 M at z ~ 0. We first explore satellite infall histories: group preprocessing 
and ejected orbits are critical aspects of satellite evolution, and properly accounting 
for these, satellite infall typically occurred at z <~ 0.5, or ~ 5 Gyr ago. To obtain 
accurate initial conditions for the SFRs of satellites at their time of first infall, we 
construct an empirical parametrization for the evolution of central galaxy SFRs and 
quiescent fractions. With this, we constrain the importance and efficiency of satellite 
quenching as a function of satellite and host halo mass, finding that satellite quenching 
is the dominant process for building up all quiescent galaxies at M star < 10 10 M Q . 
We then constrain satellite star formation histories, finding a 'delayed-then-rapid' 
quenching scenario: satellite SFRs evolve unaffected for 2 — 4 Gyr after infall, after 
which star formation quenches rapidly, with an c-folding time of < 0.8 Gyr. These 
quenching timescales arc shorter for more massive satellites but do not depend on 
host halo mass: the observed increase in satellite quiescent fraction with halo mass 
arises simply because of satellites quenching in a lower mass group prior to infall (group 
preprocessing), which is responsible for up to half of quenched satellites in massive 
clusters. Because of the long time delay before quenching starts, satellites experience 
significant stellar mass growth after infall, nearly identical to central galaxies. This 
fact provides key physical insight into the subhalo abundance matching method. 

Key words: mcthods:numerical - galaxies: clusters: general - galaxies: evolution - 
galaxies: groups: general - galaxies: haloes - galaxies: star formation. 



1 INTRODUCTION 

Observations have long shown that galaxies in denser regions 
are more likely to have low star formation rates (SFR), lie on 
the red sequence, and exhibit more evolved (elliptical) mor- 
phologies than similar mass galaxies in less dense regions, 
from massive galaxies in clusters (Oemler 1974; Davis & 
Geller 1976; Dressier 1980; Dressier & Gunn 1983; Postman 
& Geller 1984; Balogh et al. 1997; Poggianti et al. 1999) to 
the lowest mass satellites in the Local Group (Mateo 1998). 

Large-scale galaxy surveys, such as the Sloan Digital 
Sky Survey (SDSS; York et al. 2000), have enabled de- 
tailed examinations of the correlations between these galaxy 
properties and their environment at z ~ (see Blanton & 



Moustakas 2009, for a recent review). Several such early 
works showed that galaxy SFR/color depends on small- 
scale (< 1 Mpc) environment, with little-to-no additional 
dependence on larger scale environment (Hogg et al. 2004; 
Kauffmann et al. 2004; Blanton et al. 2005). More physi- 
cally, this environmental dependence has been shown to re- 
sult from satellite galaxies and the properties of their host 
dark matter halo (Weinmann et al. 2006; Blanton & Berlind 
2007; Wilman et al. 2010; Tinker et al. 2011), where 'satel- 
lite' galaxies are all those that are not the massive 'central' 
galaxy at the core of the host halo. 

These results are physically meaningful, given that the 
virial radius corresponds to a physical transition from the 
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low-density 'field' environment to a high-density, virialized 
region. After a satellite falls into a host halo, the strong grav- 
itational tidal forces prevent the satellite's (sub)halo from 
accreting dark matter and also strip mass from the sub- 
halo from the outside-in (e.g., Dekel et al. 2003; Diemand 
et al. 2007). Additionally, if the host halo is massive enough 
to host a stable, virial accretion shock (Dekel &: Birnboim 
2006), then its thermalized gas also can heat and strip any 
extended gas in the subhalo (Balogh et al. 2000; Kawata & 
Mulchaey 2008; McCarthy et al. 2008). Therefore, satellites 
(eventually) experience reduced gas cooling/accretion rates 
onto their disc after infall, a phenomenon known as 'strangu- 
lation' or 'starvation' (Larson et al. 1980). More drastically, 
in the extreme case of both high gas density and satellite 
velocity, ram-pressure can strip cold gas directly from the 
disc (Gunn & Gott 1972; Abadi et al. 1999; Chung et al. 
2009). The dense collection of galaxies in a host halo also 
allows for the possibility of strong gravitational interactions 
with neighboring galaxies, known as 'harassment' (Farouki 
& Shapiro 1981; Moore et al. 1998), and satellites can merge 
with one another (Makino & Hut 1997; Angulo et al. 2009; 
Wetzel et al. 2009b,a; White et al. 2010; Cohn 2012). All of 
these mechanisms are expected to play some role in quench- 
ing satellite star formation, though their importance, par- 
ticularly as a function of host halo mass, remains in debate. 

To constrain satellite quenching processes, many 
works have examined the SFRs/colors of satellites in 
groups/clusters in detail at z w (e.g., Balogh et al. 2000; 
Ellingson et al. 2001; De Propris et al. 2004; Weinmann et al. 
2006; Blanton & Berlind 2007; Kimm et al. 2009; van den 
Bosch et al. 2008; Hansen et al. 2009; Kimm et al. 2009; 
Pasquali et al. 2009; von der Linden et al. 2010; Prescott 
et al. 2011; Peng et al. 2012; Wetzel et al. 2012; Woo et al. 
2013). In general, these works found that the fraction of 
satellite that are quiescent/red, /q , depends primarily, and 
independently, on three quantities: /g at (1) increases with 
satellite mass, (2) increases with the mass of the host halo, 
and (3) increases toward halo center. Trend (1) is caused, 
at least partially, by the underlying dependence on stellar 
mass set by central galaxies prior to infall. Trend (2) is some- 
times interpreted as satellites being quenched more rapidly 
in more massive host halos, but the hierarchical nature of 
halo growth, namely, the possibility of quenching as a satel- 
lite in a lower mass halo prior falling into a more massive 
halo ('group preprocessing') complicates this interpretation 
(e.g., Zabludoff & Mulchaey 1998; McGee et al. 2009). Fi- 
nally, trend (3) implies an evolutionary trend, because a 
satellite's halo-centric radius negatively correlates with its 
time since infall (e.g., Gao et al. 2004; De Lucia et al. 2012). 
Similar satellite trends persists out to at least z ~ 1 (e.g., 
Cucciati et al. 2006; Cooper et al. 2007; Gerke et al. 2007; 
Tran et al. 2009; Peng et al. 2010; McGee et al. 2011; George 
et al. 2011; Muzzin et al. 2012). 

Several works have gone beyond simple satellite 
SFR/color cuts to examine observationally the nature of 
the full SFR/color distribution. These works have shown 
that the color (Balogh et al. 2004; Skibba 2009) and SFR 
(Balogh et al. 2004; McGee et al. 2011; Peng et al. 2012; 
Wetzel et al. 2012; Woo et al. 2013; Wijesinghe et al. 2012) 
distribution of galaxies is strongly bimodal across all envi- 
ronmental/halo regimes, and the SFR/color of actively star- 
forming/blue galaxies does not vary with any environmental 



measure. As noted in many of the above works, these results 
imply that the environmental process(es) takes considerable 
time (several Gyrs) to affect satellite SFR. 

In this paper, we seek to use the aforementioned obser- 
vational trends, which we presented in detail in Wetzel et al. 
(2012), to quantify — in a robust, statistical manner — the 
star formation histories and quenching timescales of satel- 
lite galaxies at z — across a wide range of both satel- 
lite and host halo masses. Understanding satellite quench- 
ing mechanisms and the timescales over which they oper- 
ate is important for elucidating the physical processes that 
occur in groups and clusters, but also for a comprehen- 
sive understanding of galaxy evolution overall, given that 
satellites constitute ~ 1/3 of all low-mass galaxies (e.g., 
Yang et al. 2007). Satellite galaxies also provide unique 
laboratories for examining gas depletion and its relation to 
star formation because, unlike central galaxies, satellites are 
thought not to accrete gas from the field after infall. Fur- 
thermore, because satellites are significantly more likely to 
lie on the red sequence, many methods for identifying galaxy 
groups/clusters rely on selecting red-sequence galaxies (e.g., 
Gladders & Yee 2000; Koester et al. 2007), so a detailed un- 
derstanding of the systematics of these methods requires 
characterizing the timescale over which satellites migrate 
onto the red sequence after infall and how this timescale 
depend on host halo mass and redshift. 

Many works have investigated satellite SFR evolution 
and quenching through the use of semi-analytic models 
(SAMs) applied to cosmological Af-body simulations. In one 
early example, Balogh et al. (2000) modeled satellite SFR as 
declining exponentially after infall on a cold gas consump- 
tion timescale of a few Gyrs to account for radial gradients 
of average SFR and color in clusters. Many SAMs assumed 
that a satellite subhalo's hot gas is stripped instantaneously 
as it passes within the host halo's virial shock, but this sce- 
nario quenches star formation too rapidly; only models that 
remove/strip gas more gradually produce realistic quiescent 
fractions (Weinmann et al. 2006; Font et al. 2008; Kang & 
van den Bosch 2008; Book & Benson 2010), but, in general 
they have difficulty in correctly reproducing the full SFR 
distribution. Though, Weinmann et al. (2010) recently im- 
plemented a modification of the SAM of De Lucia & Blaizot 
(2007) in which the diffuse gas around a satellite galaxy is 
stripped at the same rate as its host dark matter subhalo 
(10 - 20% loss per Gyr), showing that this modification pro- 
duces a satellite SSFR distribution that broadly is in agree- 
ment with observations. However, understanding the results 
of SAMs is complicated by the fact that they do not com- 
pletely accurately model the evolution of central galaxies, 
so they do not provide fully accurate initial conditions at 
infall for satellites. Relatedly, a few works have examined 
satellite SFR evolution in cosmological hydrodynamic simu- 
lations of galaxy groups (e.g., Feldmann et al. 2011), arguing 
that quenching occurs largely through the lack of gas accre- 
tion and (to a lesser degree) gas stripping after infall. 

Instead of attempting fully to forward-model all of the 
relevant physical processes for satellites, our approach is to 
parametrize satellite star formation histories and constrain 
their quenching timescales in as much of an empirical man- 
ner as possible. We start with detailed measurements of the 
SFRs of satellites at z = from our SDSS group catalog that 
we presented in Wetzel et al. (2012) and review in §2. We 
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also describe our cosmological ./V-body simulation, which we 
use to create a mock group catalog to compare our models to 
observations robustly. With this simulation, we explore the 
infall times of satellites in §3. We then develop an accurate, 
empirical parametrization for the initial SFRs of satellites at 
their time of infall in §4.1. Having accurate initial SFRs of 
satellites and measurements of their final SFRs at z — 0, we 
examine the importance and efficiency of satellite quench- 
ing in §4.2 and their star formation histories and quenching 
timescales in §4.3. With this, we then examine where satel- 
lites were when they quenched in §5.1 and their stellar mass 
growth in §5.2. Finally, we discuss the implications of our 
results for subhalo abundance matching in Appendix A. 

This paper represents the third in a series of four. In 
Tinker et al. (2011), hereafter Paper I, we described our 
SDSS galaxy sample, presented our method for identifying 
galaxy groups/clusters, and showed that central and satel- 
lite galaxy quiescent fractions are essentially independent of 
the large-scale environment beyond their host halo. In Wet- 
zel et al. (2012), hereafter Paper II, we used our SDSS group 
catalog to examine in detail the SFR distribution of satellite 
galaxies and its dependence on stellar mass, host halo mass, 
and halo-centric radius, finding that the SFR distribution is 
strongly bimodal in all regimes. Based on this, we argued 
that satellite star formation must evolve in the same man- 
ner as central galaxies for several Gyrs after infall, but that 
once satellite quenching starts, it occurs rapidly. In Wetzel 
et al. (2013a), hereafter Paper IV, we will examine quench- 
ing in galaxies near groups and clusters, focusing on 'ejected' 
satellites that passed within a more massive host halo but 
have orbited beyond the virial radius. We will show that 
these ejected satellites can explain essentially all trends for 
star formation quenching in galaxies beyond the the virial 
radius of groups/clusters. Finally, in Wetzel et al. (2013b), 
hereafter Paper V, we will use the detailed orbital histories 
from our simulation to constrain the physical mechanisms 
responsible for satellite quenching. 

For clarity, we outline some nomenclature. We refer to 
galaxies as 'quiescent' in an observational sense: having low 
SFR but without regard to how or when SFR faded. By 
contrast, we refer to satellite 'quenching' in our models as 
the physical process of SFR fading rapidly below the qui- 
escence threshold, under the ansatz that once a satellite is 
quenched it remains so indefinitely. Our galaxy group cata- 
log refers to 'group' in a general sense, as a set of galaxies 
that occupy a single host halo, regardless of its mass, and 
we will use '(host) halo' as a more general term for group 
or cluster. Finally, we cite all masses using h — 0.7 for the 
Hubble parameter. 



2 METHODS 

In this section, we first briefly describe our galaxy sample 
and group-finding algorithm. (For full details, see Paper I 
for our galaxy sample and group-finding algorithm, and Pa- 
per II for our SFR measurements.) We then describe our 
simulation, subhalo finding and tracking, and methodology 
for making galaxy and group catalogs in the simulation. 



2.1 SDSS Galaxy Catalog 

Our galaxy sample is based on the NYU Value-Added 
Galaxy Catalog (Blanton et al. 2005) from SDSS Data Re- 
lease 7 (Abazajian et al. 2009). Galaxy stellar masses are 
from the kcorrect code of Blanton & Roweis (2007), as- 
suming a Chabrier (2003) initial mass function (IMF). We 
construct two volume-limited samples of all galaxies with 
M r -5 log(ft) < -18 and -19 1 , which go out to z = 0.04 and 
0.06, from which we identify stellar mass completeness limits 
of 5 x 10 9 and 1.3 x 10 10 M®, respectively. Combining these 
samples leads to an overall median redshift of z = 0.045, 
though we will indicate this as z = for brevity. 

For our galaxy star formation metric we use specific 
star formation rate, SSFR = SFR/M sta r, based on the cur- 
rent release 2 of the spectral reductions of Brinchmann et al. 
(2004), with updated prescriptions for active galactic nu- 
clei (AGN) contamination and fiber aperture bias correc- 
tions following Salim et al. (2007). These SSFRs are de- 
rived primarily from emission lines (mostly Ha), but in 
cases of strong AGN contamination or no measurable emis- 
sion lines, the SSFRs are inferred from D n 4000. Roughly, 
SSFR > 10 -11 yr _1 are based almost entirely on Ha, 
10~ 12 < SSFR < 10" 11 yr 1 are based on a combination 
of emission lines, and SSFR < 10~ 12 yr -1 are based almost 
entirely on D„4000 and should be considered upper limits to 
the true value (Salim et al. 2007). The use of spectroscopi- 
cally derived SSFRs is critical for our analysis because dust 
reddening causes simple red/blue color cuts to overestimate 
the quiescent fraction by up to 50%, particularly at lower 
mass (see Fig. 1 in Paper I). 

2.2 SDSS Group Catalog 

Motivated by the paradigm that all galaxies reside in host 
dark matter halos, we identify groups of galaxies that oc- 
cupy the same host halo and their halo properties through a 
modified implementation of the group-finding algorithm of 
Yang et al. (2005, 2007). For our group catalog, we define 
dark matter host halos such that the mean matter density 
interior to the virial radius is 200 times the mean back- 
ground matter density: M 2 oo m = 200p m fvrTiloom- We place 
galaxies into groups through an iterative procedure outlined 
in Paper I, using the M star > 10 9 ' 7 M Q sample at z < 0.04 
and the M star > 10 10 ' 1 M sample at 0.04 < z < 0.06. 
We assign dark matter halo masses to groups by matching 
the abundance of halos above a given dark matter mass to 
the abundance of groups above a given total stellar mass: 
n(> M viri haio) = n(> M s ta r , group). Here, we use the host 
halo mass function from Tinker et al. (2008), based on a 
flat, ACDM cosmology of fi m = 0.27, fi b = 0.045, h = 0.7, 
n s = 0.95 and er 8 = 0.82, consistent with a wide array 
of observations (e.g., Komatsu et al. 2011, and references 
therein). Every group contains one 'central' galaxy, which 
by definition is the most massive, and can contain any num- 
ber (including zero) of less massive 'satellite' galaxies. 

We define a group's center by the location of its most 
massive galaxy. However, in reality the most massive galaxy 

1 In Paper II, we wrote this as M r < — 18 and -19, but we did 
not assume a value for h in calculating magnitudes. 

2 http : //www . mpa-gar ching . mpg . de/SDSS/DR7/ 
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is not always the one closest to the minimum of a halo's 
potential well, particularly at high halo mass (Skibba et al. 
2011). This effect arises largely because of the shallowness 
of the M a tar — M( sub ) halo relation at high mass, such that 
any non-negligible scatter in this relation leads to a signifi- 
cant probability that if a less massive halo falls into a more 
massive halo, the central galaxy of the less massive halo 
has higher stellar mass than that of the more massive halo. 
However, these are typically cases in which two galaxies in a 
halo have similarly high mass, and in this regime almost all 
galaxies exhibit quiescent SSFRs regardless of central versus 
satellite demarcation. Also, because we will assume realis- 
tic scatter (0.15 dex) in the M star — M( sub ) halo relation in 
making our simulation group catalog (see §2.4), this effect 
propagates into our model results as well. 

Because of projection effects and redshift-space distor- 
tions, our groups inevitably contain interloping galaxies: 
some central galaxies are mis-assigned as satellites in higher 
mass halos (reducing group purity), and conversely some 
satellites are mis-assigned as central galaxies of lower mass 
halos (reducing group completeness). As detailed in Paper I, 
an average of ~ 10% of galaxies are mis-assigned in this way 
(see also Yang et al. 2007). In this work, we apply the same 
group-finding algorithm to our simulation, as described be- 
low, allowing us to examine our theoretical results in 'obser- 
vational' space and correct for these effects. 

2.3 Simulation and Subhalo Tracking 

Our goal is to understand the SFR evolution of galaxies in 
groups/clusters and how this evolution connects with satel- 
lite infall times and lifetimes, which are governed by complex 
dynamical processes. To this end, we require a cosmological 
simulation that both provides significant statistics across a 
broad range of host halo masses and can robustly track satel- 
lite evolution from first infall to final merging/disruption. 
We employ a dissipationless, iV-body simulation using the 
TreePM code of White (2002) with flat, ACDM cosmology 
of fl m = 0.274, fit = 0.0457, h = 0.7, n = 0.95 and a 8 = 0.8, 
nearly identical to the cosmology used in making our group 
catalog. To achieve both high resolution and significant vol- 
ume, the simulation evolves 2048 3 particles in a 250 hT x Mpc 
box, with a particle mass of 1.98 x 10 s M and a Plum- 
mer equivalent smoothing of 2.5/i~ 1 kpc. Initial conditions 
are generated at z — 150 using second-order Lagrangian 
Perturbation Theory, with a displacement RMS of 38% of 
the mean inter-particle spacing. 45 outputs are stored from 
z = 10 to 0, spaced evenly in ln(a), with output time spac- 
ings of 400 and 650 Myr at z = 1 and 0, respectively. This 
same simulation was used in White et al. (2010). 

We identify 'host halos' using the Friends-of-Friends 
(FoF) algorithm (Davis et al. 1985) with a linking length 
of b = 0.168 times the mean inter-particle spacing, which 
groups particles bounded by an isodensity contour of ~ 
lOOx the mean matter density, (b — 0.2 is often used, but 
it is more susceptible to joining together distinct, unbound 
structures.) Note that this halo definition is different from 
the spherical overdensity definition used in making the SDSS 
group catalog, but we address this issue in §2.5. 

Within host halos, we identify 'subhalos' as overden- 
sities in phase space through a 6-dimensional FoF algo- 
rithm (F0F6D), also described in White et al. (2010). Based 



on extensive experimentation, we use a configuration space 
linking length of 0.078 of the simulation's mean interpar- 
ticle spacing and a velocity linking length of 0.368 of each 
host halo's ID velocity dispersion. 3 Our tests show that our 
F0F6D implementation leads to good agreement with SUB- 
FIND (Springel et al. 2001) for massive, well-resolved sub- 
halos, but F0F6D is significantly more robust in tracking 
low-mass subhalos and those that pass close to halo center. 
For both host halos and subhalos, we keep all objects with 
at least 50 particles, and we define its center and velocity 
by the position and velocity of its most bound particle. 

We track host halos and subhalos across simulation out- 
puts and build merger trees as described in Wetzel et al. 
(2009a) and Wetzel & White (2010), with slight modifica- 
tions as given below. We assign to each (sub)halo a unique 
'child' (sub)halo at the next simulation output, based on its 
20 most bound particles. We track subhalo histories across 
four consecutive outputs at a time because a subhalo can 
briefly disappear while passing through the center of its 
host halo or another subhalo. In these cases, we ensure that 
the subhalo is identified at each output by creating 'virtual' 
subhalos via interpolating the properties of the temporarily 
disappearing subhalos between outputs. If a (sub)halo has 
multiple 'parent' (sub)halos at the previous output, we iden- 
tify the main parent as the most massive one (using M max 
for subhalos, see below), and we use this main parent in 
tracking back a (sub)halo's history and identifying its main 
progenitor at an earlier time. 

We define a 'central' subhalo as being the most mas- 
sive subhalo in a newly-formed (has no parent) host halo at 
a given simulation output. The central subhalo almost al- 
ways corresponds to the object at the minimum of a halo's 
potential well. A subhalo retains its 'central' definition un- 
til falling into (specifically, becoming linked via FoF to) a 
more massive host halo, becoming a 'satellite' subhalo. Ev- 
ery sufficiently bound halo hosts one central subhalo at its 
core and can host zero, one, or multiple satellite subhalos, so 
these 'central' and 'satellite' definitions for subhalos closely 
reflect those of galaxies in the group catalog. 

We assign to each subhalo a maximum mass, M max , 
motivated by the strong correlation of this quantity with 
galaxy stellar mass (see §2.5). A subhalo's M max is based on 
the maximum host halo mass that it ever had as a central 
subhalo (so it corresponds to FoF halo mass and not F0F6D 
subhalo mass). For a central subhalo, M max almost always 
corresponds to its current halo mass, the primary exception 
being those that have passed through a more massive halo 
and been ejected (see §3.1). For a satellite subhalo, M max 
almost always corresponds to its halo mass sometime prior 
to infall, though not necessarily immediately prior to infall 
because a satellite typically undergoes some mass stripping 
just prior to infall, arising from the strong tidal forces near 
a massive halo. (The two effects above motivate our use of 
Afmax instead of mass at infall, which was used in Wetzel 
et al. (2009a,b) and Wetzel & White (2010).) Thus, a satel- 



In rare cases, using these parameters leads to zero subhalos in 
a low-mass halo. To avoid having halos with no subhalos, mostly 
for tracking purposes, we slowly increase the linking lengths in 
those halos, and we stop if we identify at least two subhalos. This 
procedure does not affect halo mass ranges used in this work. 
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lite's M max remains fixed after infall, unless it merges with 
another satellite, in which case the resultant satellite is given 
the sum of its parents' Af max values. 

As demonstrated in Wetzel &: White (2010), simulations 
at our high resolution scale can resolve and track massive 
satellite subhalos past the point at which the galactic stel- 
lar component that they host would (start to) be stripped, 
merge with the central galaxy, or otherwise be disrupted. 
Not accounting for this galactic merging/disruption in N- 
body simulations can lead to stronger small-scale clustering 
than seen in observations. This effect is even more signifi- 
cant for our phase-space F0F6D subhalo finder, which tracks 
subhalos more robustly down to small halo radii. To ac- 
count for this effect, we use the subhalo merging/disruption 
scheme in Wetzel & White (2010) and remove subhalos 
with Mbound/M max < 0.007, which provides good agree- 
ment with mass/luminosity-dependent clustering, satellite 
fractions, and luminosity functions in clusters for this sim- 
ulation and F0F6D subhalo finder. Thus, we properly re- 
solve the orbital histories and infall times of all satellites in 
our sample (M sta r > 5 x 10 9 M Q corresponds to subhalo 
M max > 3 x 10 11 M Q , see below). 



2.4 Simulation galaxy catalog 

Under the assumption that a galaxy resides at the center of 
each surviving dark matter subhalo, we assign stellar mass 
using subhalo abundance matching (SHAM; Vale & Ostriker 
2006; Conroy et al. 2006). This method assumes a one-to- 
one mapping that preserves rank ordering between subhalo 
M max (or maximum circular velocity) and galaxy M star , such 
that n(> M max ) = n(> M star ), allowing one to assign M star 
to subhalos empirically using an observed stellar mass func- 
tion (SMF) that is recovered, by design. SHAM has suc- 
ceeded in reproducing many observed galaxy statistics, in- 
cluding spatial clustering, satellite fractions, cluster lumi- 
nosity functions, and luminosity-velocity relations (Conroy 
et al. 2006; Vale & Ostriker 2006; Berrier et al. 2006; Wang 
et al. 2006; Yang et al. 2009; Wetzel & White 2010; Trujillo- 
Gomez et al. 2011). We note that, despite these successes, 
SHAM in its simplest incarnation, using M max for both 
satellite and central subhalos, may not be fully accurate in 
assigning stellar mass to both satellite and central subhalos 
simultaneously, because there is some freedom in allowing 
satellites to follow a different relation (Neistein et al. 2011; 
Yang et al. 2012; Moster et al. 2013). We discuss this issue 
further in Appendix A. 

For this work, we use the SMF from Li & White (2009), 
based on the same SDSS NYU-VAGC sample as our galaxy 
catalog, including the same 7f -correction and IMF. We ap- 
ply SHAM at a simulation output of z — 0.05, close to the 
median redshift of our SDSS catalog. While SHAM in its 
simplest implementation assumes a one-to-one correspon- 
dence between M max and M star , a scatter of 0.15 — 0.2 dex 
in this relation is suggested by observations (e.g., Zheng 
et al. 2007; Yang et al. 2008; More et al. 2009; Wetzel & 
White 2010; Leauthaud et al. 2012). Thus, in our imple- 
mentation we assume 0.15 dex log-normal scatter in M star 
at fixed M max , achieved by deconvolving the observed SMF 
with a log-normal filter such that we recover the observed 
SMF after adding this scatter. 



2.5 Simulation group catalog 

To make robust comparisons with our SDSS group catalog, 
we produce a 'simulation group catalog' by applying the 
same group-finding algorithm that we use in SDSS to our 
simulation galaxy catalog. (We use the distant observer ap- 
proximation and do not produce a light-cone.) While we base 
our models of SFR evolution on true satellite versus central 
demarcation in the simulation, we effectively 'observe' the 
results at z = 0.05 through the simulation group catalog, 
which includes the effects of interloping galaxies caused by 
redshift-space distortions and any other systematics of the 
group-finding algorithm. In Paper V, we will show that the 
galaxy distributions in the simulation group catalog closely 
match those of the SDSS group catalogs. 

Our simulation group catalog also allows us to correct 
for the effects from interloping galaxies in measuring satel- 
lite and central galaxy quiescent fractions in the SDSS group 
catalog. Interlopers have little effect on central galaxy quies- 
cent fractions because central galaxies strongly outnumber 
satellites, but interlopers do cause the observed satellite qui- 
escent fractions to be ~ 10% too low (see Appendix C of Pa- 
per I). To correct for these effects, we create a 'mock' SDSS 
group catalog by empirically assigning SSFRs to galaxies in 
the simulation group catalog, matching the observed SSFR 
distribution separately for satellite and central galaxies in 
narrow bins of galaxy and host halo mass. We then mea- 
sure quiescent fractions in the mock catalog according to 
each galaxy's true (real-space) satellite/central designation, 
which provides values unbiased by redshift-space distortions. 

Finally, our use of a simulation group catalog mitigates 
any inconsistency between the simulation's FoF halo defini- 
tion, which allows arbitrary morphology, and the spherical 
overdensity halo definition applied to SDSS, because we al- 
ways compare the two galaxy catalogs using the same halo 
definition. For reference, M2oo m ~ 1.2Mf f(6 = 0.168). 



3 SATELLITE INFALL TYPES AND TIMES 

To inform our models for the evolution of satellite SFR and 
interpret our results, we first use the simulation to explore 
the pathways and times of infall for satellites at z = 0.05 
(the median redshift of our SDSS catalog). To highlight 
physical trends free from the ambiguities of redshift-space 
distortions, in this subsection we do not use the simulation 
group catalog, but we examine satellites based directly on 
the simulation halo catalog. 

3.1 Satellite infall and ejection 

In examining the infall times of satellites, we consider two 
ways to define infall: the time of most recent infall into the 
current host halo, or the time of first infall into any host 
halo. The latter includes any time spent in a lower mass 
halo before falling into the current host halo. The latter 
also naturally incorporates satellites that temporarily orbit 
beyond their host halo's virial radius, 7? v ir, which we call 
'ejected satellites', and then fall in again (Gill et al. 2005; 
Ludlow et al. 2009; Wang et al. 2009). If one considers only 
their most recent infall, these ejected satellites might appear 
to be central galaxies falling in directly from the field for 
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Figure 1. Infall type fractions for satellites at z = 0.05 in the 
simulation versus host halo mass. Infall types are defined at the 
time of most recent infall into their current host halo. Solid red 
curve shows satellites that fell in as a central galaxy, for the first 
time, directly from the field. Short-dashed blue curve shows satel- 
lites that fell in as a satellite in a lower mass host halo (group 
infall) . Long-dashed green curve shows satellites that expirienced 
re-infall after being ejected beyond iivir- Shaded widths show beta 
distribution 68% confidence intervals. Satellite residing in higher 
mass host halos are more likely to have fallen in as a satellite or 
experienced an ejected phase. The dependence on satellite stellar 
mass is much weaker (not shown, see text). Fewer than half of 
satellites in halos > 10 14 Mq fell in directly from the field. 

the first time, but as we will show, it is more appropriate to 
consider them as satellites since their time of first infall. 

We define the 'first' infall of a subhalo as the first time 
that it becomes a satellite (is linked by the FoF algorithm) 
in another host halo, and 'recent' infall as the time that it 
falls into the main progenitor of its current host halo. In 
defining first infall, we additionally require that a subhalo 
remains a satellite for at least two consecutive simulation 
outputs to avoid cases of temporary bridging as a subhalo 
briefly passes through the outskirts of a host halo. 

Fig. 1 shows the fraction of infall types for satellites in 
the simulation at z = 0.05, as a function of and current host 
halo mass. Shaded widths show 68% confidence interval as 
given by a beta distribution (Cameron 2011). In low-mass 
halos, 90% of satellites in our mass range fell in as central 
galaxies directly from the field for the first time, but this 
fraction drops below 50% in halos > 10 14 Mq. Thus, within 
the most massive host halos, the dominant mode of infall is 
as a satellite in a lower mass halo (group preprocessing). 
So, group preprocessing of star formation is potentially im- 
portant in this regime, as we will explore in §5.1. We also ex- 
amine the dependence on satellite stellar mass (not shown), 
though it is much weaker than for host halo mass, with the 
fraction that fell in as a satellite in another halo or for the 
second (or more) time after being ejected both dropping 
from 20% at M star = 10 9 ' 7 Mq to 10% at Af star = 10 113 Mq. 
These infall fractions agree broadly with those of previous 
works (Berrier et al. 2009; McGee et al. 2009; De Lucia et al. 
2012), though each of those works find quantitative differ- 
ences arising largely from different satellite and group mass 
limits (see discussion in De Lucia et al. 2012). 



Fig. 1 also shows that satellites that have been ejected 
beyond their host halo's R v i T play an important role in satel- 
lite evolution. We find that these satellites spend typically 
1 — 2 Gyr inside a host halo when they first fall in, experienc- 
ing a single pericentric passage before being ejected. They 
then spend longer time (2 — 3 Gyr) orbiting as a central 
galaxy beyond 7? v ir until they fall in again. Both timescales 
have no significant dependence on satellite mass, but both 
do have broad distributions extending out to 4 — 5 Gyr, im- 
plying that some satellites are kicked out via multi-body en- 
counters after several orbital passages (Ludlow et al. 2009). 

After being ejected, > 90% of these (now central) galax- 
ies continue to lose halo mass, with a typical ejected satellite 
currently having half of the halo mass that it had at the time 
of ejection. 4 This halo mass stripping occurs as ejected satel- 
lites orbit in the hot, dynamic environment surrounding a 
massive host halo. This continued halo mass loss strongly 
suggests that ejected satellite galaxies evolve in a similar 
way as those within _R v i r and would also exhibit truncated 
SFRs. In Paper IV, we will examine the trends of ejected 
satellites as a function of halo-centric distance, showing that 
the ejected fraction around massive halos can account for 
the enhanced quiescent fraction of central galaxies out to 
~ 2.5i?200m around massive clusters that we noted in Pa- 
per II (see also Wang et al. 2009). Thus, we conclude that 
SFRs in ejected satellite continue evolve in a similar manner 
as those within J? v ir, and we will treat the two populations 
identically in our models for SFR evolution. 

3.2 Satellite infall times 

With the above definitions of first and recent infall, Fig. 2 
shows how surviving satellites' time since infall (and redshift 
of infall) depends on their current host halo mass. Time since 
recent infall does not depend on either host halo mass or 
satellite stellar mass (latter not shown). There is a slight 
decrease at high halo mass, because more massive halos 
formed more recently, gaining many satellites via the recent 
infall of groups (Fig. 1). By contrast, time since first infall 
exhibits a strong increase with halo mass, a result of the 
increased fraction of infalling groups and ejected satellites 
from Fig. 1 (see also Wang et al. 2007; De Lucia et al. 2012). 
We also find that more massive satellites experienced their 
first infall slightly more recently, with median time since in- 
fall falling from 5 Gyr at M star = 10 97 M to 3.5 Gyr at 
M star = 10 11 ' 3 M . This stellar mass dependence is caused 
by hierarchical structure growth, such that later infalling 
satellites had more time to grow in mass before infall, cou- 
pled with shorter lifetimes for more massive satellites, such 
that satellites that fell in early but have survived to the 
present are preferentially of lower mass. 

Half of satellites in our mass range first fell in at z > 0.5, 
with a broad tail out to z > 1 (shaded region), so they typ- 
ically have experienced > 4 Gyr evolving as a satellite. Fur- 
thermore, while a satellite typically has spent ~ 3 Gyr in 
its current host halo, the timescale is twice as long if in a 
> 10 14 M Q halo. Thus, on average, satellites have spent 1/3 - 

4 If a subhalo grows in M max by > 50% since ejection, we define 
it to be a 'newly' formed halo and discard it from the ejected pop- 
ulation, though this affects only a few percent of ejected satellites. 
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Figure 2. Median time since infall (or redshift of infall) versus 
host halo mass for all satellites in the simulation at z = 0.05. 
Blue curve shows first infall into any host halo and red curve 
shows most recent infall into the current host halo. Shaded region 
shows 68% distribution interval, which is similar in magnitude 
for first and recent infall (latter not shown). Time since recent 
infall exhibits no significant mass dependence. Time since first 
infall exhibits a strong increase with host halo mass, driven by 
the increasing fraction of infalling groups and ejected satellites in 
more massive halos (Fig. 1). Time since first infall also exhibits a 
weak decrease with stellar mass (not shown, see text). 



1/2 of their galactic lifetime evolving as a satellite, highlight- 
ing the importance of satellite evolution as a component of 
galaxy evolution. These infall times also highlight the impor- 
tance of obtaining accurate satellite SFR initial conditions 
prior to infall to understand their subsequent SFR evolution. 

Finally, we emphasize that we limit our analysis to sub- 
halos that our simulation resolves well (at least 1500 par- 
ticles at infall), and that we have accounted for satellite 
merging/disruption in a way that matches various observed 
galaxy statistics, providing accurate satellite lifetimes. Thus, 
we expect our results to be robust against resolution effects. 



3.3 Defining satellite infall 

We have shown that more than half of satellites in host ha- 
los > 10 14 M fell in as a satellite in another halo or as 
an ejected satellite, and that, as a result, time since infall is 
higher in more massive halos if considering first infall. Com- 
bined with three key observational results from Paper II, this 
provides strong evidence that the satellite-specific quench- 
ing process(es) begins at first infall, regardless of the host 
halo's mass. First, satellites exhibit an enhanced quiescent 
fraction even in the lowest mass host halos that we probed 
(3 x 10 11 Mq). To the extent that this remains true at 
higher redshift, then the satellite-specific process(es) begins 
upon infall into any halo, regardless its mass. Second, central 
galaxies out to ~ 2 i?200m around massive host halos exhibit 
an enhanced quiescent fraction, implying that ejected satel- 
lites can become/remain quenched. Third, the satellite qui- 
escent fraction increases with host halos mass even at fixed 
projected distance from halo center, d plo j/ Rvii, including 



satellites at d„ 



i? v i r , which tend to have fallen into their 



current halo recently. More specifically, at fixed d pro j/R v 



time since recent infall does not increase with halo mass, but 
time since first infall does (see Paper V). 

Taken together, these constitute strong evidence that 
satellite-specific quenching process(es) starts to set in at first 
infall, and we will use only this definition of infall henceforth. 
For a satellite that experienced first infall at time tint, we 
refer to its time since first infall, t s inccinf = t— tint- 



4 THE QUENCHING OF STAR FORMATION 
IN SATELLITES 

In this section, we present our main results on satel- 
lite quenching. First, in §4.1, we construct an empirical 
parametrization for the evolution of SFRs of central galax- 
ies to provide accurate initial conditions for the SFRs of 
satellites at their time of first infall. Then, in §4.2 we ex- 
plore the importance and efficiency of satellite quenching 
in an empirical manner by combining satellite infall times 
with our parametrization for initial quiescent fractions from 
§4.1.1. Finally, in §4.3 we develop models for satellite SFR 
evolution after infall to constrain the timescales over which 
satellites are quenched. 



4.1 SFR in satellites at the time of first infall 

To understand the evolution of SFR in satellites after infall, 
we first need accurate initial conditions for their SFRs just 
prior to infall, as given by the SFRs of central galaxies of 
the appropriate stellar mass at a satellite's redshift of infall. 
As we showed in §3.2, satellites in our mass range first fell 
in typically at z w 0.5, with a broad distribution out to 
z > 1, so using z ~ central galaxy properties as the initial 
conditions for satellites, which has been a common approach 
(e.g., van den Bosch et al. 2008; Tinker & Wetzel 2010; Peng 
et al. 2012; De Lucia et al. 2012), is a poor approximation. 

Thus, we require a statistical parametrization of the full 
SFR distribution of central galaxies as a function of both 
stellar mass and redshift. In order to describe this quan- 
tity both accurately and in a manner that is independent of 
(and minimally degenerate with) our models of satellite SFR 
evolution after infall, we proceed in an empirical manner. 
While some previous approaches have modeled the evolu- 
tion of galaxy stellar mass and SFR empirically (e.g., Conroy 
& Wechsler 2009), these studies examined only the average 
SFRs of galaxies, not taking into account the bimodal na- 
ture of the SFR distribution, in particular, quiescent galax- 
ies. Our approach is more comprehensive, as our evolution 
parametrization contains two key components: the fraction 
of central galaxies that are quiescent and the normalization 
of galaxy SFRs. We describe these in turn. 



4-1.1 Evolution of the quiescent fraction for central 
galaxies 

Our goal is to start with observations of the evolution of 
the quiescent fraction for all galaxies out to z = 1 and 
decompose this into the evolution for satellite and central 
galaxies separately, thus allowing us to use the values for 
central galaxies to provide the initial quiescent fractions for 
satellites at infall. In order to quantify reasonable systematic 
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modeling uncertainty, we implement two contrasting param- 
eterizations for this satellite-central decomposition at z > 0. 
Our 'fiducial' parametrization separates satellite and cen- 
tral galaxies based on analysis of spatial clustering measure- 
ments (Tinker & Wetzel 2010): here, the satellite quiescent 
fraction does not evolve at fixed stellar mass, and the cen- 
tral galaxy quiescent fraction declines rapidly with redshift. 
As an 'alternative' parametrization, we allow both satellite 
and central galaxy quiescent fractions to evolve at the same 
rate, which leads to the central galaxy quiescent fraction 
declining more gradually with redshift. We include both of 
these contrasting parametrizations with the goal of bracket- 
ing reasonable modeling uncertainty in our approach. 

We begin by describing our fiducial parametrization. We 
can obtain the fraction of central galaxies that are quiescent, 
fq n = nQ U /n C en, by knowing the fraction of all galaxies that 
are quiescent, fg 1 — Uq /n , the fraction of all galaxies that 
are satellites, / sat = n sat /n a11 , and the fraction of satellites 

fsat 



/sat 

that are quiescent, /q" 



reel 

JQ 



■- n^f/nsat, vi 

rail rsat r 

JQ ~ JQ Jsi 



(i) 



1 — /sat 

in which each fraction depends on stellar mass and redshift. 

To determine the quiescent fraction for all galaxies, 
/q U ( A/ S tar, z), we combine our SDSS results for all galaxies 
at z = with quiescent fractions from the COSMOS survey 
at z < 1 (Drory et al. 2009). While the results of Drory 
et al. are based on photometric spectral energy distribu- 
tions (SED), the significant number of photometric bands 
(30) in COSMOS helps to ensure accurate redshifts, stellar 
masses, and active versus quiescent demarcations. In partic- 
ular, Drory et al. identified active versus quiescent galaxies 
using full SED fitting, which minimizes the effects of dust 
contamination as compared with simpler color cuts. 5 Fig. 3a 
shows the evolution of the quiescent fraction in bins of stellar 
mass, along with the best-fit relation in each mass bin 

fg(M star , z) = A(M star ) x (1 + . (2) 

We split our SDSS sample into narrow redshift bins, and 
we anchor the fit to the quiescent fraction in the lowest 
(z < 0.04) bin. Note that the evolution within SDSS broadly 
agrees with the fits to much higher redshift. In all mass bins, 
the quiescent fraction has at least doubled since z — 1. 

To determine the satellite fraction, / sa t(M s tar, z), we 
use the simulation directly, motivated by the agreement 
(within observational uncertainty) of this quantity between 
the simulation and our SDSS group catalog at z = (see 
also Wetzel & White 2010). We use SHAM to assign stel- 
lar mass to subhalos at each redshift, using the SMF from 
SDSS (Li & White 2009) at z = 0.1 and from COSMOS 
(Drory et al. 2009) at z = 0.3,0.5,0.7,0.9, assuming 0.15 
dex Mstar — M max scatter in all cases. Fig. 3b shows the evo- 
lution of the satellite fraction in stellar mass bins. We find 
that linear growth with redshift at fixed mass, as given by 

/sat (AW , z) = B (AW ) + Bl ( M star ) 2 , (3) 

5 Despite using the same IMF, the SED-based stellar masses in 
Drory et al. are estimated to be ~ 0.2 dex higher on average than 
those in our SDSS catalog. However, their stellar masses also have 
larger scatter because of their photometric redshifts. In terms of 
quiescent fractions, these effects largely cancel out, so we do not 
attempt to renormalize stellar masses. See Tinker et al., in prep.. 
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Figure 3. (a): Fraction of all galaxies that are quiescent versus 
redshift, in bins of stellar mass. Points at z fti show our results 
from SDSS, while points at higher redshift are from the Cosmic 
Evolution Survey (COSMOS) (Drory et al. 2009). Curves show 
best-fit evolution in each mass bin. (b): Fraction of galaxies that 
are satellites versus redshift from the simulation, using SHAM to 
assign stellar mass at each redshift. Curves show best-fit evolution 
in each mass bin. (c): Fraction of central galaxies that are quies- 
cent versus redshift. Solid curves show lower limits from equation 
(1) and dashed curves show upper limits from equation (5). The 
difference from (a) is small at high mass but more significant at 
low mass. Note that the quiescent fractions in (a) and (c) mono- 
tonically increase with stellar mass, while the satellite fractions 
in (b) monotonically decrease with stellar mass. 
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Fiducial parametrization 

fcen( M /^ l (Ar S tar,*)-/^ t (M 8 ta, 

/ * J" til" n r ^ * J i x' / 71 

'Q V ' ' l-/ B at(M Btar 


)/ B at(M at ar,z) 


Fit 


Parameter 


Value 






log(Af sta r/ Mq) 




9.5 - 10.0 


10.0-10.5 10.5-11.0 11.0-11.5 



f$\z) = A{i + zy 

/sat (z) = Bq+BjZ 

/?f(M star ) = C + Ci log(M 8tar ) 



A 0.227 0.471 0.775 0.957 

a -2.1 -2.2 -2.0 -1.3 

B 0.33 0.30 0.25 0.17 

Bi -0.055 -0.073 -0.11 -0.10 

C -3.26 

Ci 0.38 



Alternate parametrization 
/§ en (M 8 ta r , z) = /§ en (M star , z = 0) x (1 + zY 



f^(M stal ,z = 0) = Do + D 1 log(M Bta . T ) D 0(): , 



Do -6.04 
i 



Table 1. Fits of the quiescent fraction for central galaxies, as a function of stellar mass and redshift. Top: Fiducial parametrization, 
based on equation (1), which incorporates the quiescent fraction for all galaxies, fq (z), the satellite fraction, / sa t(.z), and the quiescent 
fraction for satellites, /Q at (M s t ar ). Fits to /q (z) and /sat(-z) are given in stellar mass bins, while /^ at (Af a tar) is fit to stellar mass 
dependence assuming no redshift dependence. Bottom: Alternate parametrization, based on equation (5). /g en (M s tar) is fit to stellar 
mass dependence at z = 0, and a, in bins of stellar mass, is the same as in the fiducial parametrization. 




Figure 4. Fraction of galaxies that are (were) quiescent versus 
stellar mass. Solid curves show values at z = measured in the 
SDSS group catalog (using the simulation group catalog to correct 
for interloping galaxies), for satellite (blue, top) and central (red, 
middle) galaxies. Dotted lines show fits, given in Table 1. Green 
region at bottom shows what fraction of satellites at z = were 
quiescent prior to their first infall, with uncertainty given by our 
two parametrizations. The significant increase of satellites' qui- 
escent fraction from first infall (green) to z = (blue) highlights 
the importance of satellite quenching. The difference between the 
quiescent fraction for satellites prior to first infall (green) and for 
central galaxies at z = (red) shows that the latter does not 
represent accurate initial conditions for satellites. 



provides a reasonable fit, as Fig. 3b shows. Note that the 
fluctuations with redshift are driven by the evolution of the 
SMF in Drory et al. and not by subhalo statistics. 

Finally, to determine the quiescent fraction for satel- 



lites, 



we first fit to the stellar mass depen- 



dence at z = from our SDSS group catalog, as Fig. 4 
shows. Recall that we use the simulation group catalog to 
remove the effects of galaxy interlopers caused by redshift- 
space distortions. We find that 

/q *(Mrt«, z = 0) = Co + Ci log(M star ) (4) 

provides a reasonable fit in our mass range, as shown by the 
dotted curve. We then impose that there is no evolution of 
this quantity, such that /^(AW, z) = /cf (M star ,2 = 0). 
This choice is motivated by the results of Tinker & Wet- 
zel (2010), who found no evolution in the quiescent frac- 
tion for satellites at fixed magnitude at z < 1, based on 
halo occupation modeling of the spatial clustering and num- 
ber densities of galaxy samples from the Classifying Objects 
by Medium-Band Observations (COMBO-17) (Phleps et al. 
2006) and Deep Extragalactic Evolutionary Probe (DEEP2) 
(Coil et al. 2008) surveys. Note that Tinker et al., in prep, 
find similar results from the spatial clustering of COSMOS 
galaxies. 

Putting these ingredients into equation (1), we obtain 
the quiescent fraction for central galaxies as a function of 
stellar mass and redshift, as Fig. 3c shows (solid curves). 
The difference from the overall quiescent fraction (panel a) 
is modest at high mass but is more significant at low mass, 
where central versus satellites quiescent fractions differ more 
strongly and the satellite fraction is higher. 

This fiducial parametrization is based on the satellite 
quiescent fraction not evolving at fixed stellar mass since 
z = 1. This behavior is motivated by spatial clustering mea- 
surements and also is supported by George et al. (2011), 
who examined galaxy groups of mass 10 13 " 14 M in COS- 
MOS out to z = 1 and found no significant evolution in the 
quiescent fraction of group members, at least for sufficiently 
massive galaxies (A/ S t ar > 3 x 10 10 Mq) that their sam- 
ple is complete. However, several other works observe that 
the quiescent/red fraction of galaxies in clusters decreases 
with increasing redshift (e.g., Butcher & Oemler 1984; Pog- 
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gianti et al. 2006; McGee et al. 2011). In order to investigate 
this possible systematic uncertainty, we develop an alternate 
parametrization to quantify the impact of non-zero evolution 
of the satellite quiescent fraction at fixed stellar mass. Mo- 
tivated by the idea that, at least out toz~ 1, the quiescent 
fraction for satellites is always higher than for central galax- 
ies of the same mass, as supported by observations (e.g., 
Cooper et al. 2007; McGee et al. 2011; George et al. 2011), 
we parametrize the quiescent fractions for both satellite and 
central galaxies as decreasing with redshift at the same rate 
as for all galaxies, such that /Q en //Q at remains fixed. 6 That 
is, at fixed mass, central galaxies evolve as 

f Q an (M st ^, z) = f Q cl \M s ^, z = 0) x (1 + z) a(M ^ , (5) 

using the same value of a as for all galaxies in the same 
stellar mass bin. We fit for /Q en (M star , z = 0) directly from 
the group catalog, as Fig. 4 shows, finding that 

fq n (M star , z = 0) = D + Di log(M star / M ) (6) 

provides a reasonable fit. The dashed curves in Fig. 3c 
show the resultant evolution of the central galaxy quies- 
cent fraction in this alternate parametrization, which ex- 
hibits a more gradual decline with redshift. However, for 
both parametrizations, the central galaxy quiescent fraction 
increases by at least a factor of two from z = 1 to 0. 

Again, we emphasize that our two parametrizations pro- 
vide significant contrast, and we consider their difference as 
representative of the reasonable systematic uncertainty in 
satellite initial conditions. 7 Table 1 lists the fits and pa- 
rameters for each term in equations (1) - (6). For binned 
parameters, we will use spline interpolation as a function of 
log (Mstar) to obtain smooth stellar mass dependence. 

4-1.2 Evolution of SFR for central galaxies 

Having parametrized the evolution of the quiescent fraction 
for central galaxies, we now develop a prescription for the 
evolution of their SFRs, which are observed to increase with 
redshift (e.g., Noeske et al. 2007). For central galaxies that 
remain active at z = 0, we parametrize their star formation 
history through a modified exponential r model, 

SFR ccn (t) oc (t-^expj-^^H (7) 

M star (i) = /retain / SFR(t) dt 
Jt f 

with tf being the time of initial formation, which we take to 
be t(z — 3) for all galaxies in our mass range, and /retain be- 
ing the fraction of stellar mass that is not lost through super- 

6 In detail, /Q° n //Q at must evolve somewhat because / sa t evolves, 
but /sat evolution causes /Q Cn //Q at to decrease by < 10% to z = 1. 

7 We also tried an extreme scenario in which satellites are solely 
responsible for the evolution of the quiescent fraction for all galax- 
ies, to see if it is possible that the quiescent fraction for central 
galaxies does not evolve at fixed mass. However, even in the ex- 
treme scenario of no quiescent satellites by z = 1, the central 
galaxy quiescent fraction at fixed mass still must decrease by a 
factor of at least 2 from z = to 1 to account for the overall galaxy 
quiescent fraction — there are simply not enough satellites. 



novae and stellar winds, which we take to be /retain = 0.6. 8 
To obtain r con , we place all active central galaxies in our 
SDSS group catalog into narrow (0.2 dex) bins of stellar 
mass and compute the median r CO n in each bin using equa- 
tion (7) with the stellar mass and median SSFR of the bin. 
This yields T CO n values that range from 3.8 to 1.9 Gyr from 
Mstar = 5 x 10 9 to 2 x 10 11 M Q . We use the median r C c n 
of each bin to evolve back the SFRs of all active galaxies in 
that bin, which agree with measured SFRs at 0.2 < z < 1.1 
from Noeske et al. (2007) to within 0.15 dex across our mass 
range, well within their measurement errors. This prescrip- 
tion keeps the width of the active galaxy SFR distribution 
fixed, also in agreement with Noeske et al.. Thus, by design, 
our r model agrees well with the full active galaxy SFR 
distributions from Noeske et al. out to z = 1. While our r 
model is highly simplified compared with the actual star for- 
mation histories of galaxies, we note that similar r models 
have been shown to agree with observed average SFRs out 
to z = 1 (e.g., Noeske et al. 2007b). Thus, for our specific 
purpose of assigning statistically accurate SFRs to active 
satellites at their time of infall, we consider this simplified 
but constrained model a reasonable empirical approach. 

Parametrizing the possible change in the SFRs of qui- 
escent central galaxies is more difficult, given both noisier 
measurements at z = and a lack of detailed SFR measure- 
ments at higher redshift. However, these uncertainties are 
largely irrelevant for us in practice, because satellites that 
fall in being already quiescent are at or quickly evolve be- 
low SSFR w 10" 12 yr" 1 , where our SDSS measurements are 
largely upper limits. For completeness, in our parametriza- 
tion we assume that quiescent galaxies at all redshifts have 
the same SFR normalization as quiescent galaxies at z = 0. 
We also explored letting the SFR normalization for quies- 
cent galaxies evolve by the same amount as for active galax- 
ies, such that the separation between the peaks in the SFR 
distribution remains unchanged, though doing this does not 
alter significantly our model results at z = in §4.3.2. 

Note that equation (7) ignores the contribution from 
galaxy mergers when computing stellar mass growth, which 
means that it underestimates the amount of total stellar 
mass growth since z — 1 somewhat. We examined the im- 
portance of this effect by using the stellar masses of galax- 
ies from SHAM out to z = 1 and following galaxy-galaxy 
merger histories in our simulation to z = 0. The typical 
amount of stellar mass growth via mergers since z — 1 for 
galaxies in our stellar mass range is always < 10% (usually 
much less); mergers become important only for galaxies of 
significantly higher mass. Thus, this introduces only a small 
bias in our parametrization as compared with modeling un- 
certainty of initial quiescent fractions in §4.1.1. 

In summary, to obtain the SFR distribution of central 
galaxies as a function of mass and redshift we first compute 
the median T ccn for active central galaxies at a given stel- 
lar mass, we use this r con to compute the increase in SFR 
to a given redshift, and we apply this increase to all active 



In our model, this mass loss occurs instantaneously, though in 
reality, it is an extended process. As Leitner & Kravtsov (2011) 
showed, the vast majority (~ 90%) of stellar mass loss occurs 
within the first 2 Gyr, so our approximation is good given that 
90% of satellites fell in earlier than 2 Gyr ago. 
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Figure 5. Example of the evolution of the SSFR distribution 
for central galaxies at fixed stellar mass. The relative fraction of 
active and quiescent galaxies at each redshift is given by Table 1 
(fiducial parametrization) , and the SSFR normalization of active 
galaxies increases with redshift according equation (7), using the 
median r ccn of active central galaxies in the mass range. 



central galaxies at that mass, such that the active galaxy 
SFR distribution width does not evolve at fixed mass. We 
also assume that the SFR normalization for quiescent cen- 
tral galaxies does not evolve. We then evolve the relative 
fraction of active and quiescent central galaxies according 
to equation (1) or (5) by moving randomly chosen SFR val- 
ues from the quiescent to the active side of the distribution. 

Fig. 5 shows an example of the resultant evolution of 
the SSFR distribution, for central galaxies with M star — 
10 io-io.s Mq We emphasize that 

our methodology repro- 
duces the observed evolution of the quiescent fraction as 
well as the active galaxy SFR normalization and distribu- 
tion width, providing accurate (statistical) initial conditions 
for the SFRs of satellites prior to first infall. We also note 
that, while our parametrization clearly also has interesting 
implications for the physics of how central galaxies evolve, 
we defer further investigation in this area to future work. 
For this paper, our parametrization is important only in 
providing accurate initial conditions for satellite SFRs. 



4-1-3 Assigning SFR to satellites at mfall 

We now describe how we assign initial SFRs to satellites at 
their time of first infall. For each satellite in the simulation 
at z — 0, we have its time of first infall, but in order to use 
our above parametrization, we must also know what stellar 
mass each satellite at z = had at that time, at least in a 
statistically average sense. To estimate this, we posit that 
satellites have grown in stellar mass by the same amount, on 
average, as active central galaxies of the same stellar mass. 
This ansatz allows us to use equation (7), with the same me- 
dian value of T ccn as for central galaxies of the same stellar 
mass at z = 0, to calculate the average factor by which satel- 
lites were less massive at their time of first infall. As we will 
demonstrate in §5.2, this ansatz is not only self-consistent, 
but moreover it must be satisfied in our parametrization, at 



least in the absence of significant, systematic stellar mass 
loss from tidal stripping. 9 

Note that this approach, based on statistical aver- 
aged star formation histories, ignores scatter in stellar mass 
growth. As a check on the accuracy and self-consistency of 
our approach, we also tried computing stellar mass growth 
for central galaxies directly from SHAM, by differencing 
each central galaxy's stellar mass at z = from what it 
had at a given redshift according to SHAM. While this al- 
ternate approach leads to higher scatter in the amount of 
stellar mass growth over a given time interval (for example, 
20—30% since z = 1), it leads to average stellar mass growths 
that are consistent with the above method to within 10%. 
Thus, our approach is self-consistent, at least in parametriz- 
ing the average stellar mass growth of active galaxies. 

Thus, for each satellite (including ejected satellites) in 
the simulation at z = 0, we assign its SFR at its redshift of 
first infall by drawing randomly from the central galaxy SFR 
distribution at the appropriate stellar mass at that redshift, 
using equation (7). For satellites that fell in prior to z = 1, 
we extrapolate the quiescent fraction and r model fits to 
higher redshift, though most satellites that fell in at z > 1 
are low-mass and had minimal likelihoods of being quiescent 
at infall, so changing the extrapolation changes our results 
by only a few percent. 

To highlight the importance of our parametrization for 
assigning accurate initial quiescent fractions to satellites, the 
green region in Fig. 4 shows the fraction of satellites at z — 
that were quiescent prior to first infall. The region bound- 
aries are determined by our two parametrizations, with solid 
and dashed curves corresponding to those in Fig 3c. If con- 
volved with satellite infall times, the resultant quiescent- 
prior-to-infall fractions differ by less than 10%, small com- 
pared with the marked difference from the quiescent fraction 
for central galaxies at z — (red curve), which has been 
assumed for satellite initial conditions in previous works, as 
mentioned above. Furthermore, the significant difference be- 
tween satellite quiescent fractions at infall and at z = (blue 
curve) indicates the importance of the satellite quenching 
process, which we will explore next. 



4.2 Importance & efficiency of satellite quenching 

We now explore the importance of satellite star formation 
quenching in building up the full population of quiescent 
(red-sequence) galaxies at z — as well as the efficiency by 
which satellites are quenched. The results in this subsection 
are essentially empirical, relying only on our parametrization 
for satellite initial quiescent fractions from §4.1.1. These re- 
sults are independent of any particular model for the mech- 
anism^) or timescale of satellite quenching but provide in- 
sight into the efficiency of the mechanism(s) as a function 
of stellar mass. 



9 If surviving satellites have experienced significant, systematic 
stellar mass loss from tidal stripping, this would mean that they 
were more massive at infall, so their initial quiescent fractions 
were higher, than in our model. However, for feasible amounts of 
average stripping of < 30% (see Appendix A), Fig. 4 (green re- 
gion) shows that the initial quiescent fractions would not increase 
by more than ~ 5%. 
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Figure 6. Importance of satellite quenching in building up the 
quiescent population at z = as a function of stellar mass. 
Top, red region shows the fraction of currently quiescent satel- 
lites that quenched as satellites, "Qassat/ n Q t now i Middle, blue 
region shows the fraction of all currently quiescent galaxies that 
quenched as satellites, "Q t assat / tl Q 1 now > Region widths indicate 
uncertainty in satellite initial quiescent fractions from §4.1.1. Be- 
low M s tar = 10 10 Mq, satellite quenching is the dominant route 
for building up the entire quiescent (red-sequence) population. 
Dotted blue curve shows the (incorrect) n??' i/n?! 1 from 

v ' Q as sat' now 

using central galaxies at z = for satellite initial conditions. 

The basic quantity that we use in this subsection is the 
number density of satellites (including ejected satellites) at 
z — that quenched as satellites, given by 

"■qL sat (M star (Z = 0)) = riQ* now (M st ar (z = 0)) - (8) 

^inf (Mstar (^ = 0)) . 

Here, nQ* now is the number density of satellites that are qui- 
escent at z — from the SDSS group catalog (again, using 
the simulation group catalog to correct for interloping galax- 
ies), and fiQjnf is the number density of satellites at z — 
that were quiescent prior to infall from our initial condition 
parametrization. We express the latter quantity as a func- 
tion of stellar mass at z — 0, based on our assumption that 
satellites have grown in stellar mass by the same amount as 
central galaxies (which we will justify in §5.2). 

4-2.1 Importance of satellite quenching 

We first examine the contribution of satellite quenching 
to building up the quiescent population at z = 0, as 
Fig. 6 shows. First, to understand the importance of satel- 
lite quenching on just the satellite population, the red re- 
gion shows what fraction of currently quiescent satellites 
were quenched as satellites, iiQ t aS sat/ n Q'now The region 
width indicates the uncertainty from our two initial con- 
dition parametrizations in §4.1.1. At M s tar < 10 10 Mq, es- 
sentially all quiescent satellites quenched as satellites. At 
Mstar ~ 10 11 M©, this fraction is half, because half of satel- 
lites were already quiescent as central galaxies prior to infall 
(Fig. 4). Thus, in our mass regime satellite quenching always 
dominates the production of quiescent satellites. 

To demonstrate the importance of satellite quenching 



on the entire galaxy population, the blue region shows what 
fraction of all currently quiescent galaxies quenched as satel- 
lites, WQljssat/ n Q n now This fraction decreases significantly 
with stellar mass, such that the vast majority of galaxy 
quenching occurs via central galaxies at M s tar > 10 11 M©. 
In this regime, central galaxies are as likely as satellites to 
be quiescent at z — (Fig. 4) , and they outnumber satellites 
by a factor of > 6 (Fig. 3b), so central galaxies dominate the 
production of the quiescent population. By contrast, at the 
low mass end, even though central galaxies still outnumber 
satellites by a factor of ~ 3, satellite quenching is so much 
stronger that satellites dominate the production of the quies- 
cent population. Thus, at M a tar < 10 10 Mq, satellite quench- 
ing is the dominant route for the build-up of all quiescent 
(red-sequence) galaxies. Furthermore, no central ('isolated') 
galaxies are observed to be quiescent at M s tar < 10 9 M© 
(Geha et al. 2012), so satellite quenching is the only process 
for quenching galaxies at such low mass. 

We emphasize the importance of accurate satellite ini- 
tial conditions for these results. Several previous works have 
attempted to infer the impact of satellite quenching by as- 
suming that satellite initial conditions can be approximated 
via central galaxies at z — (e.g., van den Bosch et al. 
2008; Tinker & Wetzel 2010; Peng et al. 2012; De Lucia 
et al. 2012). However, this assumption necessarily under- 
estimates the importance of satellite quenching, because 
central galaxies were less likely to be quenched at higher 
redshift (Fig. 3c). To highlight the impact of this assump- 
tion, the dotted blue curve in Fig. 6 shows the resultant 
n Qassat/ n Q U now if one assumes (incorrectly) that central 
galaxy SFRs at z = represent satellite initial conditions, 
that is, using 7iQ* lnf = /Q° I r iow n sat in equation (8). This as- 
sumption underestimates the true importance of satellite 
quenching by a factor of at least 50% in our mass range. 

4-2.2 Efficiency of satellite quenching 

We next examine the efficiency by which satellites are 
quenched, as given by the fraction of satellites that were 
active at infall (able to be quenched) that then quenched 
as satellites after infall, rig assat/^Ainf ■ Fig- 7 (green re- 
gion) shows this fraction as a function of stellar mass, with 
the width again indicating the uncertainty in satellite initial 
quiescent fractions. At low mass (Mstar < 10 10 M©) half of 
satellites that were active at the time of infall have been 
quenched by now, while the other half still actively form 
stars. By contrast, at high mass (M s tar > 10 11 M©) essen- 
tially all initially active satellites have been quenched. Thus, 
more massive satellites are quenched more efficiently. Phys- 
ically, this implies that more massive satellites are quenched 
more rapidly, as we will show in §4.3.1. 

To elucidate the differing quenching efficiencies for 
satellite versus central galaxies, the orange region in Fig. 7 
shows what fraction of active-at-infall satellites would have 
quenched had they instead remained central galaxies. This 
fraction is given by (/§ ml now - /§f inf ) //I'now, with /§f now 
and /q^ow being the fractions of satellite and central galax- 
ies, respectively, that are quiescent at z = 0, and /™„ 
1 — /^ I rlow . Like satellites, central galaxies also quench more 
efficiently at higher mass, though with a lower overall effi- 
ciency. From Fig. 7, it might naively appear that the differ- 
ing quenching efficiency for satellite versus central galaxies 
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Figure 7. Satellite quenching efficiency versus stellar mass. Top, 
green region shows the fraction of active-at-infall satellites that 
were quenched as satellites after infall, ^ s S at/ n A *inf ■ R e gi° n 
widths indicate uncertainty in satellite initial quiescent fractions 
from §4.1.1. More massive satellites are quenched more efficiently. 
For reference, orange region shows fraction of active-at-infall 
satellites that would have quenched had they instead remained 
central galaxies. Dotted curve shows the satellite quiescent frac- 
tion excess, f^ XCCBB , given by equation (9). 

is stronger at lower mass, though one must interpret these 
fractions carefully. In Paper II, we argued that a robust com- 
parison is given by the satellite quiescent fraction excess 

rsat rccn 

/•sat 7Q,now /Q, now /q\ 

JQ excess Teen ' \ ) 

•1 A, now 

which represents the excess fraction of satellites that were 
quenched after infall that would not have been quenched 
had they remained central galaxies. As shown by the dot- 
ted curve in Fig. 7, /Q a CXCC ss is independent of stellar 
mass. Furthermore, in Paper II we showed that the stellar 
mass independence of /^excess persists across all host halo 
masses and halo-centric radii. 

If the same physical mechanism(s) that quenches central 
galaxies also operates on satellites, then fo"l xcess indicates 
how much more of an effect the satellite-specific quenching 
mechanism(s) has. In this scenario, the invariance of fo'l xoeaB 
suggests that the efficiency of the satellite quenching pro- 
cesses) is independent of stellar mass. However, the physical 
processes that are thought to quench central galaxies — such 
as virial shock heating, mergers, active galactic nuclei — and 
their relative importance as a function of stellar mass remain 
topics of active investigation. Thus, it is unclear if such cen- 
tral galaxy quenching processes are important for satellites, 
and if they are, whether they occur before or after the on- 
set of any satellite-specific processes. We will investigate the 
physical mechanisms of satellite quenching in more detail in 
Paper V, and we note that the results in this paper do not 
depend on the exact mechanism (s) at play. 

To summarize this empirically-motivated subsection: 

(1) satellite quenching dominates the production of quies- 
cent satellites at all masses we probe, and it dominates the 
production of all quiescent galaxies at M s tar < 10 10 Mq, and 

(2) more massive satellites are quenched more efficiently. 



4.3 SFR evolution of satellites 

We now use the satellite infall times from our simulation to 
extend the results of the previous subsection and constrain 
satellite SFR evolution and quenching timescales. This sub- 
section presents the primary results of this paper. 

The relative importance of various mechanisms for 
quenching satellites — such as strangulation, ram-pressure 
stripping, and harassment — and the details of how their ef- 
fects propagate to influencing satellite star formation remain 
topics of active investigation. Nonetheless, for most plausi- 
ble physical processes, the likelihood that a satellite has been 
quenched increases with its time since infall. 

Motivated by this idea, we proceed under the follow- 
ing ansatz: if a satellite was active at the time of first in- 
fall, the parameter that determines if it has been quenched 
is simply its time since first infall, tsmceinf, the time that it 
has spent as a satellite. As we will show in Paper V, this 
simple ansatz yields satellite quiescent fractions that have 
the correct dependencies on both halo-centric distance and 
halo mass, as compared with our observational results in 
Paper II, because both the satellite quiescent fraction and 
isinceinf increase with decreasing halo-centric distance. This 
agreement implies that the scatter between quenching likeli- 
hood and tsinco inf must be small, because a scenario in which 
quenching likelihood and t s inccinf have large scatter would 
lead to satellite quiescent fraction radial gradients that are 
too shallow (see Paper V for more). 10 Thus, our ansatz pro- 
vides a well-motivated means to constrain statistically the 
timescales over which satellites are quenched. 

In this subsection, we examine satellite quenching 
timescales in two ways, the first being simpler and more 
empirical, the second being more physical. First, in §4.3.1, 
we consider quenching simply in the binary sense, that is, 
when SSFR falls below the 10 -11 yr _1 bimodality threshold. 
We constrain the t s ince inf at which satellites are quenched in 
this binary sense, which we refer to as the satellite quench- 
ing time, tQ. Then, in §4.3.2, we decompose this quenching 
time into two more physically informative timescales: the 
time delay after infall at which star formation starts to be 
quenched, tQ, delay, and the characteristic e-folding time over 
which SFR fades once quenching has started, tq, fade- 

4-3.1 Quenching timescale of satellites 

We first constrain the time since infall at which satellites are 
quenched (fall below SSFR = 10~ n yr _1 ), t Q . To do this, 
we select all satellites in the simulation at z — 0, and we use 
our satellite initial quiescent fraction parametrization from 
§4.1.1 to compute whether each was active or quiescent prior 
to infall. Those that were quiescent remains so thereafter. 
For those that were active at the time of infall, we desig- 
nate the ones with tsince inf > tQ as having been quenched. 
Using narrow bins of both satellite and host halo mass, we 

10 More rigorously, because we use a step-function threshold 
in tsince inf i selecting the maximally oldest surviving satellites 
to quench, any scatter in the relation between quenching and 
tsince inf would lead to a characteristic quenching time (at which 
50% of active-at-infall satellites have quenched) that is necessarily 
shorter. Thus, our quenching timescales are strictly upper limits, 
but, as we have argued, the scatter should be small. 
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Figure 8. Satellite quenching timescales versus stellar mass. 
Top: Time since first infall at which satellites are quenched (fall 
below SSFR = 10 -11 yr _1 ), tQ, in bins of their host halo's mass 
at z = 0. Region widths indicate uncertainty in satellite ini- 
tial quiescent fractions from §4.1.1. Bottom: Decomposing the 
above tQ into the delay time after infall before satellite quench- 
ing starts, delay: an d subsequent e-folding time over which 
SFR fades, tq_ f a( j c . Low-mass satellites take significantly longer 
to be quenched than those at higher mass, and their SFR fades 
more rapidly once quenching starts. However, satellite quenching 
timescales do not depend on the mass of their host halo. The 
uncertainty on tq_ f at j e is dominated by uncertainty in fitting the 
full SSFR distribution. 



adjust tQ until the satellite quiescent fraction in the simula- 
tion group catalog matches that of the SDSS group catalog. 
Repeating this procedure in each mass bin yields tQ as a 
function of both satellite mass and host halo mass. 

Fig. 8 (top) shows how tQ depends on satellites' cur- 
rent stellar mass, in bins of current host halo mass. Region 
widths indicate the uncertainty from satellite initial quies- 
cent fractions (solid and dashed curves correspond to those 
in Fig. 3c). As the quenching efficiency results of §4.2.2 sug- 
gested, more massive satellites are quenched more rapidly. 
The most massive satellites are quenched ~ 2 Gyr after in- 
fall, while those at M sta r < 10 10 Mq form stars actively for 
~ 5 Gyr after infall before being quenched. Note that 5 Gyr 
is over half the age of the Universe at the typical redshift 



that they fell in, so low-mass satellites have spent as much 
as half of their entire star-forming lifetimes as satellites. 

In Paper II, we showed that satellites at z = are more 
likely to be quiescent if they reside in more massive host 
halos. But interestingly, in Fig. 8 we find no significant, sys- 
tematic dependence of satellites' tQ on their current host 
halo mass (there are fluctuations at low stellar mass, but 
they are not monotonic). This halo mass independence nat- 
urally arises from our tying satellites' quenching to their 
time since first infall, which incorporates the increased in- 
fluence of group infall and ejection/re-infall in more massive 
halos, leading to a natural increase in f s i n coinf with halo mass 
(Fig. 2). Recall that this choice was motivated by the ab- 
sence of a minimum host halo mass at z = for affecting 
satellite SFR. To the extent that this fact holds true out to 
z ~ 0.5 (~ 5 Gyr ago, the timescale over which tQ is sensi- 
tive), Fig. 8 indicates that there is little-to-no freedom for tQ 
to depend on host halo mass, because we obtain tQ by em- 
pirically matching the observed quiescent fraction at z = 
in each host halo mass bin. In other words, given that a 
~ 10 12 Mq host halo quenches a low-mass satellite ~ 4 Gyr 
after infall, as demanded by Fig. 8, a similar halo that then 
falls into a massive cluster will bring in satellites that al- 
ready are quenched. Thus, the overall tQ in massive clusters 
is set by a combination of tQ from this group preprocessing 
and from satellites that fell directly into the cluster from 
the field, but as Fig. 8 shows, this overall tQ is effectively 
the same for massive clusters as for a Milky- Way halo. 1 
Thus, the timescale over which satellites are quenched does 
not depend on the mass of their host halo. 

This lack of dependence on host halo mass may be sur- 
prising, given that more massive host halos represent more 
severe environments, having higher gas densities, tempera- 
tures, and satellite orbital velocities at a given scaled dis- 
tance, d/R v i r . But it is not clear that all possible satellite 
quenching processes should depend on host halo mass. For 
example, if satellite quenching is driven simply by the in- 
ability to accrete gas after infall, then quenching occurs 
when a satellite exhausts its gas reservoir, independent of its 
host halo's mass. Alternately, while the dominant satellite 
quenching process(es) may be more rapid at a given d/i? v ir 
in a more massive host halo, this is mitigated by the fact 
that dynamical friction causes a satellite of a given mass to 
orbit to smaller d/R v [ T more quickly in a lower mass host 
halo (Boylan-Kolchin et al. 2008; Jiang et al. 2008). We will 
examine the dependence of satellite quenching on host halo 
mass with physically motivated models applied to orbital 
histories in Paper V. 



4-3.2 ' Delay ed-then-rapid' quenching of satellites 

While the satellite binary quenching timescale that we mea- 
sured above, tQ, is advantageous in its simplicity, it is in- 
sensitive to the details of how satellite SFR evolves. We now 



11 Given that our SDSS group catalog constrains iq a function 
of host halo mass at z = and not as a function of host halo 
mass at the time of infall directly, if group preprocessing were 
the primary mode of quenching satellites, this could mitigate the 
inferred dependence of tQ on host halo mass. However, as we will 
show in §5.1, most satellites quenched when they were in their 
current host halo, so any possible mitigation would be modest. 
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seek to understand satellite star formation histories more 
fully, as constrained by the full SSFR distribution at z = 0. 

In Paper II, we showed that the satellite SSFR distri- 
bution is bimodal, similar to central galaxies, across our 
stellar mass range. The SSFR values of the active galaxy 
peak and bimodality break, as well as the fraction of galax- 
ies near the bimodality break ('green valley'), do not vary 
with central versus satellite demarcation, host halo mass or 
halo-centric radius. As we argued, these observations imply 
that (1) satellite SFRs evolve in the same manner as cen- 
tral galaxies for several Gyrs after infall, (2) the time since 
infall at which satellite SFR starts to be affected is long com- 
pared with the time over which SFR fades, and (3) the latter 
timcscale does not depend on host halo mass or halo-centric 
radius. 

To quantify these timescales, we build on these trends 
and construct a physically motivated, two-stage model for 
satellite SFR evolution. The initial SFR for a satellite at its 
time of first infall, ti n f , is given by our parametrization in 
§4.1. If a satellite was quiescent prior to infall, we do not 
evolve its SFR. If a satellite was active at infall, we allow 
its SFR after infall to fade gradually in the same manner as 
central galaxies of the same stellar mass, using equation (7). 
This central-type, gradual fading continues across a 'delay' 
time, t Q> delay If isincoinf > £q, delay, only then does a satel- 
lite start to be quenched, at which point we parametrize its 
SFR evolution via exponential fading, with tq, fade being the 
characteristic e-folding time over which SFR fades. Defining 
£Q, start = iinf + tQ, delay, satellite SFR evolves as 

{SFR C en(t) * < t Q . start 

( C-'Q, start) ] 
SFR C e„(tQ, s tart)el T Q> fade / t >t Q . st!lrt 

(10) 

to the redshift of our group catalog. Note that tQ, delay 
and t q , f a dc relate to tQ from the previous subsec- 
tion Via tQ = tQ, delay + N TQ , fade, With N = 

In [SSFR(tQ, start )/10- n yr" 1 ]. 

Because the initial SFRs that we assign to satellites are 
based on observed distributions, any possible measurement 
uncertainty propagates into our resultant model SFRs at z = 
as well, allowing for robust comparison with SDSS. Also, 
for any satellite whose SSFR evolves below « 10~ 12 yr -1 , 
we assign it as having SSFR = 1CF 12 yr -1 plus log-normal 
scatter of 0.2 dex, which effectively mocks the measurement 
limits and scatter in the Brinchmann et al. (2004) method. 

Our physical, two-stage model for satellite SFR evo- 
lution has two free parameters (timescales) to constrain: 
tQ, delay and t q , f a dc ■ We allow these timescales to vary, in- 
dependently, with both satellite and host halo mass, again 
constrained to yield the observed satellite quiescent frac- 
tion at z — in the joint mass bins. The mere existence 
of a bimodal distribution with non-zero population at inter- 
mediate SSFRs requires that both timescales are non-zero. 
However, to explore the impact of the two timescales, we 
consider three scenarios: 

(a) delay = 0, tq, fade fit to the quiescent fraction 

(b) tq, fade = 0, tQ, delay fit to the quiescent fraction 

(c) tQ, delay, tq, fade jointly fit to full SSFR distribution 

Fig. 9 shows each resultant SSFR distribution at z = 
in bins of stellar mass, for all satellites in host halos with 



M2oom > 10 1 M©, using our fiducial parametrization for 
satellite initial quiescent fractions. (Using our alternate 
parametrization shifts tQ, delay to slightly longer values, as 
Fig. 8 (bottom) shows, but does not affect tq, f a de or the 
quality of our fits.) We discuss each scenario in turn. 

First, we examine scenario (a), in which satellite 
quenching begins immediately at infall (tQ. delay = 0), and 
SFR fades slowly over a long tq, f a de- Fig. 9a shows the resul- 
tant SSFR distributions that yield the correct satellite qui- 
escent fractions. This slow-fade scenario puts far too many 
satellites at intermediate SSFRs, violating the bimodality 
and leading to a qualitatively incorrect distribution. Satel- 
lite SFR does not quench gradually after infall. 

Next, we consider the opposite scenario (b), in which 
satellite quenching is delayed after infall, but once it starts, 
it occurs instantly (tq, f a dc = 0). Fig. 9b shows the resul- 
tant SSFR distributions, which have a qualitatively correct 
bimodality, including a correct SSFR values of the active 
peak and bimodality break. In particular, the agreement of 
the SSFR distribution for active galaxies confirms that they 
have evolved in the same manner as active central galaxies. 
However, the bimodality break is clearly too strong, with 
a deficit of galaxies at intermediate SSFRs. Satellite SFR 
must take non-trivial time to fade once quenching starts. 

Finally, in scenario (c) we allow both tQ, delay and 
TQ,fade to vary, as fit to the full SSFR distribution, pro- 
viding a unique solution, as Fig. 9c shows. This simple, two- 
stage quenching scenario produces a SSFR distribution in 
excellent agreement with observations at each stellar mass 
bin, with tq, f a de being 10 — 20% of tQ, delay In Fig. 8 (bot- 
tom), we show more explicitly how the best-fit tQ, delay plus 
tq, fade times from Fig. 9c depend on current satellite mass, 
in bins of current host halo mass, including uncertainty in 
tQ, delay from satellite initial quiescent fractions. As with tQ 
in the previous subsection, both tQ, delay and tq, f a de do not 
depend on the mass of the host halo. We do not plot sep- 
arate curves for tq, fade in bins of host halo mass, because 
the dominant uncertainty in tq, f ade , shown in Fig. 8, comes 
from fitting the full SSFR distribution, which is larger than 
any systematic changes with host halo mass or satellite ini- 
tial conditions. Furthermore, as mentioned in §4.1.2, we find 
that these tq, fade values do not depend on whether or not we 
evolve the SFR normalization for quiescent central galaxies. 

To summarize this subsection, satellite SFR evolves, at 
least on average, via a 'delayed-then-rapid' quenching sce- 
nario: satellite SFR remained unaffected for 2 — 4 Gyr af- 
ter first infall (depending on stellar mass), after which SFR 
fades rapidly, with an e-folding time of < 0.8 Gyr. Both 
timescales are shorter for more massive satellites have no 
significant dependence on host halo mass. 



5 IMPLICATIONS OF SATELLITE 
QUENCHING TIMESCALES 

Using the quenching timescales that we constrained in 
§4.3.2, we now explore two implications for satellite evo- 
lution. First, in §5.1 we explore where satellites were when 
they quenched, focusing on the importance of group prepro- 
cessing. Second, in §5.2 we use our constrained model for 
SFR evolution to examine how much satellites have grown 
in stellar mass since infall. 
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Figure 9. Specific star formation rate (SSFR) distribution of satellites at z = 0, in bins of stellar mass, in all host halos with Af200m > 
10 12 Mq. Solid blue curves show results from the SDSS group catalog while dashed red curves show results from the simulation group 
catalog using the satellite SFR evolution model of equation (10) under the following scenarios, (a): Satellite SFR starts to be quenched 
immediately upon infall and fades over an e-folding time, tq i fade- This scenario produces no bimodality. (b): Satellite SFR remains 
unaffected after infall for a quenching delay time, tQ, delay: after which it quenches instantaneously. This scenario produces a bimodality 
that is too strong, (c): Satellite SFR remains unaffected after infall for time iq, delay: after which it fades over an e-folding time TQ^fade- 
This 'delayed-then-rapid' quenching scenario produces the correct SSFR distribution. These results are insensitive to host halo mass, 
and varying satellite initial quiescent fractions affects only iQ, delay (see Fig. 8). 



5.1 Where were satellites when they quenched? 

We have argued that the physical process(es) responsible 
for quenching satellites sets in at the time of first infall, but 
that its effects take considerable time to propagate before 
quenching starts. We now ask: where were satellites at the 
moment when they started quenching? In §3.1, we examined 
what fraction of satellites fell in directly from the field versus 
as a satellite in a group. We now extend those results, using 
the quenching times from §4.3, to examine what fraction of 
currently quiescent satellite quenched (a) prior to first infall, 
(b) in a different host halo prior to falling into their current 
host halo, or (c) in their current host halo? 

For each surviving satellite in the simulation, we com- 
pute if it was active at the time of first infall, as before. If so, 
we use the tQ t delay values from §4.3.2 for the time at which 
it started to quench. We then compute whether the satellite 
was in the main progenitor of its current host halo or was 
in a different host halo at the time that quenching started. 
Fig. 10 shows what fraction of currently quiescent satellites 
quenched in the three different regimes, as a function of cur- 
rent stellar mass, in bins of current host halo mass. 

Fig. 10a shows what fraction of currently quiescent 
satellites already were quenched as a central galaxy prior 
to first infall. As discussed in §4.2, this fraction increases 
with stellar mass, both because the central galaxy quiescent 
fraction is higher at higher stellar mass and because higher 
mass satellites fell in more recently, when central galaxies 
were more likely to be quiescent. Fig. 10a also shows that 
quenching prior to first infall is more important in lower 
mass host halos, again because satellites in lower mass host 
halos fell in more recently. 

Fig. 10b shows what fraction of currently quiescent 
satellites started quenching in a different host halo prior 
to falling into their current host halo, indicating complete 
group preprocessing. Opposite to the trends for quenching 
prior to first infall, this fraction is higher for lower mass 



satellites and in higher mass host halos, both trends a result 
of the hierarchical nature of halo growth (§3). In particular, 
half of all low-mass (M sta r < 10 10 Mq) quiescent satellites 
in massive clusters (M2oom > 10 14 Mq) started quenching 
as a satellite in a group. 

Finally, Fig. 10c shows what fraction of currently quies- 
cent satellites quenched while in their current host halo. This 
mode of quenching dominates at most masses, though the 
fraction is alway < 70%; its importance wanes both at high 
stellar mass, where quenching prior to first infall dominates, 
and at low stellar mass, where the importance of group pre- 
processing increases. Only about half of quiescent satellites 
within massive clusters (> 10 14 Mq) quenched there. 

In summary, group preprocessing has a critical impact 
on satellite star formation histories. We have argued that 
any time spent as a satellite in another host halo is im- 
portant as far as starting the quenching process, but these 
results demonstrate the impact of complete group prepro- 
cessing. This is particularly important for satellites in clus- 
ters, in which 15 — 50% of all quiescent satellites started 
quenching as a satellite in another host halo. Given the hi- 
erarchical nature of halo growth, group preprocessing should 
be only more important for quenching satellites below our 
5 x 10 9 M stellar mass limit. 



5.2 Stellar mass growth after infall 

So far, we examined satellite star formation histories with a 
focus on star formation rate evolution, but in this last sub- 
section we examine the implications for stellar mass growth. 
In §4.3, we showed that satellites continue to form stars ac- 
tively, in the same manner as central galaxies, for 2 — 4 Gyr 
after infall, which represents as much as half of their total 
star-forming lifetimes. Thus, satellites have the capacity to 
grow significantly in stellar mass via star formation after 
infall. For now, we ignore any other processes that might 
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Figure 10. Where satellites were when they quenched: the fraction of currently quiescent satellites that quenched in different host halo 
regimes versus stellar mass, in bins of current host halo's mass. Region widths indicate uncertainty in satellite initial quiescent fractions 
from §4.1.1. (a): Fraction that were quenched as a central galaxy prior to first infall. Quenching prior to infall is more important in regimes 
where satelitcs fell in more recently: higher mass satellites and in lower mass host halos (§3.2). (b): Fraction that started quenching in 
a different host halo (group) prior to falling into their current host halo, based on delay from §4.3.2. Half of quiescent satellites at 
Mstar < 10 10 M that are currently in massive clusters began quenching as a satellite in a group, highlighting the importance of 'group 
preprocessing', (c): Fraction that quenched in their current host halo. Overall, quenching prior to infall dominates at high satellite mass, 
group preprocessing is significant at low satellite mass, and quenching within the current host halo dominates at intermediate mass. 
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Figure 11. Satellite stellar mass growth after infall. Top: Median 
ratio of a satellite's stellar mass at z = to what it had at 
first infall versus current stellar mass, for all satellites (middle, 
grey), those that remain active (bottom, blue) and those that are 
quiescent (top, red) at z = 0. Region widths show uncertainty in 
satellite initial quiescent fractions from §4.1.1. Bottom: Median 
ratio of a satellite's stellar mass at z = to what it would have 
at z = had it not quenched after infall. Satellite quenching has 
little impact on stellar mass growth: satellite and central galaxy 
stellar mass growths via star formation are nearly identical. 

affect satellite stellar mass evolution, such as tidal stripping 
or merging, though we discuss these in Appendix A. 

To quantify the amount of stellar mass growth via star 
formation that satellites at z = have experienced since first 
infall, we use our model for satellite SFR evolution given by 
equation (10), with the appropriate values of r ccn and the 



quenching timescales tQ, delay and tq, f a d e from §4.3.2 given 
a satellite's stellar mass at z = 0. We integrate SFR(t) to 
obtain the stellar mass formed since first infall, again assum- 
ing that 40% of this stellar mass is lost through supernovae 
and stellar winds. For a satellite that was quiescent prior 
infall, its SFR has been sufficiently low that we can neglect 
any stellar mass growth since that time. We then examine 
statistical trends by computing the median fractional stellar 
mass growth since first infall in bins of stellar mass at z — 0. 

Fig. 11 (top) shows the median ratio of a satellite's stel- 
lar mass at z = to the mass that it had at the time of its 
first infall, as a function of its current stellar mass. As be- 
fore, region widths show uncertainty in satellite initial qui- 
escent fractions from §4.1.1. Considering all surviving satel- 
lites (grey region), their median stellar mass growth since in- 
fall is negligible at high mass but is 50% at M star < 10 10 M . 
This mass dependence arises because lower mass satellites 
are more likely both to have fallen in earlier when SFRs were 
higher and to have been active at the time of infall. 

The blue and red regions in Fig. 11 show median val- 
ues for currently active and quiescent satellites, respec- 
tively. Overall, currently active satellites have experienced 
significantly less stellar mass growth since infall than cur- 
rently quiescent galaxies. While perhaps counter-intuitive, 
this trend is readily understandable. Even though active 
satellites are still growing in stellar mass, to remain active 
they necessarily fell in more recently, meaning both lower 
SFRs at the time of infall and less time for mass growth af- 
ter infall. To understand currently quiescent satellites, note 
that they are composed of two populations: those that were 
quiescent prior to infall and those that quenched after in- 
fall. While those that were quiescent prior to infall did not 
grow in stellar mass at all, those that quenched after infall 
necessarily fell in early, when SFRs were much higher, and 
they then spent several Gyrs actively forming stars before 
being quenched. While the former population dominates at 
high mass, the latter dominates at low mass, leading to stel- 
lar mass typically having more than doubled since infall at 
Mstar < 10 10 Mq. Thus, low -mass satellites experience con- 
siderable stellar mass growth after infall. 
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To put this result in context, we compare the stellar 
mass growth experienced by satellite versus central galaxies. 
Fig. 11 (bottom) shows the median ratio of satellite mass at 
z — to what it would be if all satellites that were ac- 
tive at infall remained active to z = 0, that is, if satellites 
never quench. This approximately indicates the ratio of stel- 
lar mass that satellites have compared to what they would 
have if they had remained a central galaxy. (In fact, the 
masses of satellite versus central galaxies are closer than in 
Fig. 11, because some central galaxies have quenched since 
the time that a satellite fell in, but we do not attempt to fully 
model central galaxy SFR evolution here.) While this mass 
ratio is unity for currently active satellites (by definition), 
considering the entire satellite population, the median ra- 
tio remains consistent with unity (grey region). Considering 
just currently quiescent satellites, which have experienced 
the most truncated mass growth, the median reduction in 
stellar mass at z = is never more than 10% (red region). 
Furthermore, considering just those satellites that quenched 
after infall, the reduction is still only 10%, though it remains 
at that level across all stellar mass (not shown). Thus, de- 
spite the clear importance that satellite quenching has on 
instantaneous SFR, satellite-specific quenching has minimal 
impact (< 10%) on satellite stellar mass growth; satellite 
and central galaxy stellar mass growth via star formation is 
nearly identical. This behavior arises for two reasons: satel- 
lites evolve for considerable time (2 — 4 Gyr) after infall 
until they start to be quenched, and galaxies in our mass 
range form the vast majority of their stars at high redshift 
(z > 0.5) when their SFRs were much higher, so quenching 
at low redshift has little impact on their final stellar mass. 

As outlined in §4.1, in order to assign accurate initial 
SFRs to satellites at their time of first infall, we estimated 
their stellar mass at that time via the ansatz that they 
grew by the same amount as central galaxies of the same 
stellar mass. This is not obviously a good approximation, 
for instance, if observations had constrained both delay 
and tq, fade in §4.3.2 to be quite short. However, the re- 
sults of this subsection reassuringly show that our approach 
is statistically self-consistent to good approximation. More- 
over, we investigated alternative scenarios in which satel- 
lites have grown significantly less in stellar mass since in- 
fall than central galaxies, but this generically leads to even 
longer quenching times, and thus more implied stellar mass 
growth, so these scenarios are not internally self-consistent. 
Thus, based on Fig. 11, the only self-consistent scenario is 
that stellar mass growth in satellites is the same as central 
galaxies to within 10% (at least in the absence of significant, 
systematic stellar mass loss from tidal stripping). 

Finally, the results of this section provide physical in- 
sight and possible improvement into the implementation of 
the SHAM method in assigning galaxy stellar mass to both 
central and satellites subhalos. The basic idea of SHAM is 
that one can assign instantaneous stellar mass to all subha- 
los, central and satellite, based simply on some measure of 
their subhalo mass (or circular velocity). Our results show 
that, because satellites and central galaxies grow in stellar 
mass by essentially the same amount, on average, it is jus- 
tifiable to assign stellar mass to all subhalos under a single, 
simple prescription. However, our results do imply possible 
tension with the way that SHAM typically is implemented, 
specifically, through the use of the maximum/infall subhalo 
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Figure 12. Summary diagram of satellite galaxy SFR evolu- 
tion, as given by equation (10), highlighting the 'delayed-then- 
rapid' quenching scenario. The physical process(es) responsible 
for quenching star formation in a satellite begins after it first falls 
into another host halo, regardless of the host halo's mass. How- 
ever, it takes considerable time for SFR to be affected: satellite 
SFR evolves after infall in the same manner as central galax- 
ies for a delay time t q i delay = 2 — 4 Gyr (red curve) , depend- 
ing on stellar mass. Then, the satellite's star formation starts 
to be quenched, and SFR fades rapidly, with an e-folding time 
T Q,fade = 0.2 — 0.8 Gyr, also depending on stellar mass. Less 
massive satellites have longer t q f delay 

and TQ y fadci but neither 
timescale dependends on the mass of the host halo. Because of 
the long tQ t delayi satellite stellar mass growth via star formation 
is nearly equal to that of central galaxies. Because of hierarchical 
halo growth, many satellites in massive host halos were quenched 
as a satellite in a lower mass halo prior to infall. For comparison, 
blue curve shows gradual SFR fading for central galaxies, with 
characteristic fading time T cen = 2 — 4 Gyr from equation (7). 



mass, which does not evolve after infall for satellites. We 
discuss this issue in Appendix A. 



6 SUMMARY AND DISCUSSION 
6.1 Summary 

Using a galaxy group/cluster catalog from SDSS Data Re- 
lease 7, together with a cosmological A^-body simulation to 
track satellite orbits, we examined in detail the star forma- 
tion histories of satellite galaxies at z « 0, focusing on their 
times since infall, quenching timescales, and stellar mass 
growth after infall. Applying the same group-finding algo- 
rithm to our simulation as we used in SDSS allows us to 
make robust comparisons of model results to observations. 
To obtain accurate initial conditions for the SFRs of satel- 
lites at their time of first infall, we constructed an empiri- 
cally based, statistical parametrization for the evolution of 
central galaxy SFRs out to z = 1; this is critical for the ac- 
curacy of our results because, at fixed stellar mass, the qui- 
escent fraction for central galaxies more than doubles from 
z = 1 to 0. Our primary result is that, at least on average, 
satellite SFR evolves via a 'delayed-then-rapid' quenching 
scenario: satellite SFR remained unaffected for several Gyrs 
after first infall, after which quenching occurs rapidly, as 
Fig. 12 summarizes. In more detail, our main results are: 
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Infall as part of a group and ejection beyond R v i T are 
important aspects of satellite evolution. Fewer than half of 
satellites in massive clusters fell in directly from the field; the 
rest fell in as a satellite in another host halo or experienced 
secondary infall after becoming ejected. Satellites at z — 
experienced their first infall typically at z ~ 0.5, or ~ 5 Gyr 
ago, with a broad tail out to z > 1. Less massive satellites 
and those in more massive host halos fell in earlier. 

Satellite quenching is a critical process of galaxy evolu- 
tion. Satellite quenching always dominates the production 
of quiescent satellites. Moreover, satellite quenching is re- 
sponsible for producing the majority of all quiescent (red- 
sequence) galaxies at M sta r < 10 M by z = 0. 

Satellite quenching is delay ed-then-rapid. Based on ob- 
servations, we argued that the process(es) responsible for 
quenching satellites begins after first infall into any other 
host halo. As constrained by the satellite SSFR distribution 
at z = 0, satellites, at least on average, then remain actively 
star-forming for 2 — 4 Gyr (depending on stellar mass) after 
first infall, unaffected by their host halo, before quenching 
starts. Once quenching has started, the e-folding time over 
which SFR fades is < 0.8 Gyr. 

Satellite quenching timescales are shorter at higher stel- 
lar mass but are independent of host halo mass. More mas- 
sive satellites start to be quenched more rapidly after in- 
fall, and once quenching has started, their SFRs fade more 
quickly. However, these quenching timescales do not depend 
on host halo mass. The observed increase in the satellite qui- 
escent fraction with host halo mass arises because of the in- 
creased importance of group preprocessing and ejection/re- 
infall in more massive host halos. 

Group preprocessing plays a critical role in quenching 
satellites. Half of low-mass (M atar < 10 10 M ) quiescent 
satellites in clusters started quenching in another host halo 
before falling into the cluster. Across all satellite masses, the 
fraction of quiescent satellites that were quenched within 
their current host halo is never more than ~ 70%. 

Satellite quenching barely impacts stellar mass growth. 
Because satellite quenching is so delayed, low-mass satellites 
have experienced considerable stellar mass growth via star 
formation since infall: satellites at M sta r < 10 10 M are, 
on average, 50% more massive than at infall. Moreover, the 
average amount of mass growth via star formation in satel- 
lite and central galaxies is identical to within 10%. This 
provides key physical insight into the abundance matching 
technique for assigning stellar mass to subhalos, as outlined 
in Appendix A. 

6.2 Relation to satellite gas content 

We first discuss the relation of our results to satellite gas 
content. Satellites provide unique laboratories for examin- 
ing gas depletion and its relation to star formation because, 
unlike central galaxies, satellites' subhalos are thought not 
to accrete matter after infall: the strong gravitational tidal 
forces in the host halo both prevent a satellite's subhalo 
from accreting new matter and strip any existing subhalo 
matter, including gas, from the outside-in. Additionally, any 
thermalized gas in the host halo can heat and ram-pressure 
strip any extended subhalo gas, and in the extreme case of 
both high gas density and satellite velocity, ram-pressure 
can strip cold gas directly from the disc. 



We first discuss the implications of our iQ, delay results, 
that is, that SFR in satellites evolves for 2—4 Gyr after infall, 
depending on stellar mass, unaffected by the host halo. While 
this timescale may represent, to some degree, the statistical 
average of those of individual satellites, it is informative to 
consider in the context of gas depletion times. Given that 
star formation is fueled by cold, molecular gas (Wong & 
Blitz 2002; Bigiel et al. 2008), this means that a significant 
quantity of such gas must persist in a satellite's disc for that 
amount of time. 

One possibility is that a sufficient reservoir of cold gas 
was present in the disc at the time of infall. As a constraint, 
we compare our tq. delay times to observed cold gas depletion 
times, defined as M gas /SFR. At z = 0, observed atomic gas 
depletion times in M s t ar > 10 10 M galaxies are ~ 3 Gyr, 
with large scatter but no systematic dependence on stellar 
mass or SFR (Schiminovich et al. 2010). Incorporating the 
additional ~ 30% of the gas that is molecular (Saintonge 
et al. 2011), the total gas depletion time would extend to 
~ 4 Gyr. This timescale can be even longer to the extent 
that gas recycled from stellar mass loss fuels star formation 
(e.g., Leitner & Kravtsov 2011). If valid at higher redshift, 
this depletion time would be sufficient to accommodate our 
tq, delay values if all of the atomic gas converts to stars. 

Observed gas ratios, M gas /M s t a r, provide another con- 
straint. In §5.2, we showed that satellites experience signif- 
icant stellar mass growth via star formation. In particular, 
currently quiescent, low-mass satellites have more than dou- 
bled their stellar mass since infall, which requires a gas reser- 
voir comparable in mass to their stars at the time of infall. 
Observations at z = show that actively star-forming galax- 
ies at M s tar ~ 10 10 M have total cold gas masses that are 
~ 40% of their stellar mass (Catinella et al. 2010; Saintonge 
et al. 2011), and that this gas ratio increases with decreas- 
ing stellar mass, being near unity for galaxies just below our 
mass threshold (e.g., Geha et al. 2006). However, currently 
quiescent, low-mass satellites fell in at higher redshift (typi- 
cally, z — 0.5 — 1), and if the total cold gas fraction increases 
with redshift at a rate suggested by observations of molecu- 
lar gas in actively star-forming, massive (M star > 10 10 M ) 
galaxies at z = 0.4 — 1.4 (Daddi et al. 2010; Tacconi et al. 
2010; Geach et al. 2011), then these satellites would have 
had enough cold gas in their disc to accommodate the sig- 
nificant stellar mass growth in Fig. 11. 

Furthermore, the cold gas in the disc could be replen- 
ished for some time after infall if the most concentrated and 
tightly-bound component of the extended subhalo gas con- 
tinues to cool/accrete onto the disc for several Gyrs before 
being stripped. X-ray observations show that roughly half 
of massive satellites in groups and clusters retain extended, 
hot gas halos, though truncated as compared with central 
('field') galaxies (Sun et al. 2007; Jeltema et al. 2008). Sim- 
ulations also show retention of extended subhalo gas: Mc- 
Carthy et al. (2008) found that satellite subhalos can retain a 
significant fraction (~ 30%) of their hot gas for several Gyrs 
after infall, while Simha et al. (2009) and Keres et al. (2009) 
found that satellites continue to accrete significant gas onto 
their disc, though at a reduced rate compared with central 
galaxies. Adding this replenishment to what cold gas was al- 
ready in the disc at infall, the total gas reservoir in/around 
satellites appears fully sufficient to fuel their extended SFR 
and stellar mass growth as demanded by tQ, delay 
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Finally, regarding £q, delay, we note that the halo ra- 
dius crossing time, given our virial definition, is t CIOSS — 
Rvir/Vvir — 2.7(1 + z)~ 3 ^ 2 Gyr, independent of host halo 
mass. More precisely, numerically integrating satellite orbits 
in an NFW potential (assuming energy and angular momen- 
tum conservation) using typical initial orbital parameters 
from Wetzel (2011) for satellites in host halos in our mass 
range, the average time from infall to first pericentric pas- 
sage is somewhat shorter at 2(1 + z)~ 3 ^ 2 Gyr, independent 
of satellite mass. Thus, the onset of satellite quenching oc- 
curs near the time of pericentric passage for more massive 
satellites and 1 — 2 Gyr after pericentric passage for lower 
mass satellites, as we will explore in more detail in Paper V. 

We next discuss the implications of our tq, fade re- 
sults, that is, that once quenching has started, satellite 
SFR fades rapidly, with tq, f a de being 0.8 to 0.2 Gyr from 
M st ar = 5 x 10 9 to 2 x 10 11 M©. A lack of star forma- 
tion implies a lack of dense, molecular gas, so it is interest- 
ing to compare tq, f a de to observed molecular gas depletion 
times, defined as Mh 2 /SFR, for galaxies near the quenching 
threshold. Saintonge et al. (2011) examined the molecular 
gas depletion times in a large sample of SDSS galaxies, find- 
ing typically ~ 1 Gyr at M sta r = 10 10-11 M . However, 
they found that the depletion time increases with decreas- 
ing SSFR, being ~ 2 Gyr at SSFR ~ 10 -11 yr _1 . Such a 
long depletion time in galaxies near the quenching thresh- 
old is difficult to understand in a scenario in which satel- 
lites simply use up their molecular gas, which may suggest 
that additional processes are at play in reducing the molec- 
ular gas density in satellites as they are quenched, possibly 
via tidal or ram-pressure stripping or internal feedback pro- 
cesses. However, note that the sample in Saintonge et al. 
is composed primarily of central ('field') galaxies, and it is 
unclear whether the significant gas reservoirs in these nearly 
quiescent galaxies result from them not fully depleting their 
cold gas while quenching, or by them accreting gas after 
they have quenched. If the latter holds, the cold gas proper- 
ties of satellites as they are being quenched could be quite 
different. 

Several works have examined the gas content of satel- 
lites in the Virgo cluster, finding that while the atomic gas 
masses of Virgo satellites are significantly lower than for 
galaxies of the same mass in the field (Huang et al. 2012; 
Serra et al. 2012), the molecular gas masses are quite simi- 
lar (Kenney & Young 1986; Young et al. 2011). This result 
suggest that, while processes like tidal or ram-pressure strip- 
ping may play a role in removing atomic gas from the outer 
regions of satellites, they have little impact on the molec- 
ular gas that fuels star formation. Indeed, observations of 
ram-pressure stripping (Sun et al. 2007; Chung et al. 2009; 
Abramson et al. 2011) typically show diffuse, atomic gas 
being stripped from the outer regions of the disc, while the 
dense, molecular gas towards the core survives intact, a phe- 
nomenon also seen in simulations (Tonnesen & Bryan 2009). 

Overall, we conclude that the gas reservoir in/around 
satellites at their time of infall is sufficient to fuel their nec- 
essary star formation histories and stellar mass growth. Our 
result that satellite quenching can be parametrized simply 
by time since first infall, with no significant dependence on 
host halo mass, suggests that simple gas depletion ('strangu- 
lation') most naturally explains satellite quenching, though 
more work is needed to understand if gas self-depletion 



alone can account for our 'delayed-then-rapid' quenching 
scenario with sufficiently short tq, f a de- Our tq, f a de values 
are marginally-to-significantly shorter than observed molec- 
ular gas depletion times, particularly for galaxies near the 
quenching threshold, possibly suggesting that some process 
other than simple molecular gas depletion is at play, though 
it is not clear if external stripping processes can explain this. 
In Paper V, we will examine these issues in more detail by 
developing physical models for satellite SFR evolution. 

6.3 Comparison with other work 

Our satellite quenching timescales are broadly consistent 
with previous works that parametrized the evolution of 
SFR in satellites and argued that it is affected over long 
(2 — 3 Gyr) timescales (e.g., Balogh et al. 2000; Wang et al. 
2007; McGee et al. 2009; Mahajan et al. 2011). Recently, De 
Lucia et al. (2012) examined the infall times of satellites in 
a SAM applied to the Millennium simulation, accounting for 
hierarchical halo growth and group preprocessing; compar- 
ing with observed satellite quiescent fractions as a function 
of host halo mass and halo-centric radius, they argued that 
satellites take ~ 5 — 7 Gyr to quench after falling into halos 
> 10 13 Mq. However, these previous works generally only 
used observed quiescent/red fractions to constrain a single, 
overall quenching timescale, as in our Fig. 8. A significant 
aspect of our approach is using the overall SSFR distribu- 
tion to constrain satellite SFR evolution more completely, 
through which we have shown that satellites experience a 
'delayed-then-rapid' quenching scenario, which is not pos- 
sible measuring just quiescent fractions. Furthermore, our 
use of empirically based satellite initial SFRs and a mock 
simulation group catalog to compare robustly with observa- 
tions allows us to constrain these timescales empirically and 
robustly. 

Our results place strong constraints on semi-analytic 
approaches to modeling the physics of satellite SFR evolu- 
tion (e.g., Font et al. 2008; Kang & van den Bosch 2008; 
Weinmann et al. 2010; Kimm et al. 2011). Our results sug- 
gest that a successful physical model would allow satellite 
SFR to evolve, environmentally unaffected, for 2 — 4 Gyr 
(depending on stellar mass) after infall, possibly through 
continued accretion/cooling of extended subhalo gas. Our 
results also suggest that one does not need to impose any 
explicit dependence on host halo mass to this process. 

We emphasize that our quenching timescales are valid 
for satellites at z ~ 0, given that our approach is sensitive 
to satellite SFR evolution within the last ~ 4 Gyr, and it is 
not clear that these timescales remain fixed at higher red- 
shift. We have checked this in our framework by applying our 
quenching timescales from z = to satellites at higher red- 
shift and comparing with the observed evolution in the qui- 
escent fraction for all galaxies from Fig. 3. While our model 
does agree within observational uncertainties at z < 0.3, we 
find that this approach quenches too few galaxies at higher 
z. One one hand, this discrepancy could be interpreted as 
evidence against the accuracy of our 'delayed-then-rapid' 
quenching scenario. Alternately, satellite quenching times 
simply may be shorter at higher redshift. Using halo occu- 
pation modeling of galaxy spatial clustering measurements, 
Tinker & Wetzel (2010) showed that the satellite quiescent 
fraction does not evolve with redshift at fixed magnitude. 
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This lack of satellite evolution is supported broadly by ob- 
servations of massive galaxies (M sta r > 3 x 10 10 Mq) in 
X-ray-selected groups of mass 10 13 ~ 14 M Q out to z ~ 1 m 
COSMOS (George et al. 2011). If the satellite quiescent frac- 
tion does not evolve at fixed satellite and host halo mass, 
then the ratio of a satellite's quenching time to its dynam- 
ical friction lifetime must remain roughly fixed, implying 
that the satellite quenching time is shorter at higher red- 
shift: t Q oc (1 + z)~ 3/2 (Tinker & Wetzel 2010). Further- 
more, because the SSFR distribution of galaxies in groups 
is strongly bimodal out to at least z ~ 0.4 (McGee et al. 
2011), this suggests that the delayed-then-rapid quenching 
scenario remains true at higher redshift. Also, if tg t delay de- 
creases more quickly than tq, fade, then there would be a 
higher fraction of satellites at intermediate SSFRs at higher 
redshift, a trend that is suggested by the significant fraction 
of 'green valley' galaxies in groups at z ~ 1 (Balogh et al. 
2011). In all, these result suggest that satellite quenching 
times are shorter at higher redshift, as we will investigate 
and quantify further in Paper V. 

We also emphasize that our results are predicated on the 
accuracy of the SHAM method for assigning stellar mass to 
both satellite and central subhalos across 2 = 0—1. While 
numerous works support the (statistical) accuracy of this 
approach, as outlined in §2.4, to any extent that SHAM 
assigns biased stellar masses to satellite subhalos at these 
redshifts, this would bias our derived quenching timescales. 

We have argued that group preprocessing is the primary 
reason that satellites in more massive host halos are more 
likely to be quiescent. Group preprocessing should manifest 
itself via an increased quiescent fraction for satellites that 
remain in a sub-group after falling into a cluster (White 
et al. 2010), to the extent that such sub-groups are obser- 
vationally identifiable (Cohn 2012). Promisingly, examining 
the red fraction of galaxies in clusters at z = 0.2 — 0.5, Li 
et al. (2009) saw that, within and at the outskirts of these 
clusters, galaxies that appear associated with groups exhibit 
a higher red fraction than those that are not, providing di- 
rect evidence for the importance of group preprocessing. 

Our results also connect with satellites of much lower 
mass in the Local Group. Naively extending our quenching 
timescales in Fig. 8 to lower mass implies that dwarf satel- 
lites in the Milky Way take at least 5 Gyr to quench after 
infall. Given that the vast majority of satellites in the Local 
Group are quiescent, this implies that they fell in > 5 Gyr 
ago, as supported by detailed comparisons of their positions 
and velocities to those of similar satellites in simulation 
(Rocha et al. 2012). Conversely, our quenching timescales 
also reinforce recent claims that the Large and Small Mag- 
ellanic Clouds, which are both actively star-forming, fell into 
the Milky Way halo within the last few Gyrs (Besla et al. 
2007): had they fallen in much earlier, our results indicate 
that they would no longer be actively star-forming. 

Examining the stellar age distributions of satellites in 
the Local Group, Orban et al. (2008) found that, even for 
those that are currently quiescent, a large fraction have ex- 
perienced a significant amount of recent star formation: half 
of satellites with M sta r > 10 7 Mq have formed more than 
10% of their stellar mass in the last 2 Gyr. This trend, ob- 
served in lower mass satellites than we examined, supports 
our general result that satellites continue to form stars over 
extended timescales and thus grow in stellar mass consider- 



ably after infall. Examining more massive (M sta r > 10 9 Mq), 
quiescent elliptical/lenticular galaxies in the Coma cluster, 
Trager et al. (2008) found that their mean ages are identical 
to those of similar galaxies that are in the field, and that 
they have experienced star formation as recently at z ~ 0.2, 
trends that again suggest significant star formation after in- 
fall. Finally, Smith et al. (2012) examined the stellar ages 
of quiescent satellites in the Coma cluster, finding a signifi- 
cant decrease in age with cluster-centric radius at low mass 
(Mstar < 10 10 Mq) but a much weaker trend at higher mass. 
As they argued, this trend implies that satellite-specific 
quenching plays a stronger role in quenching lower mass 
satellites, consistent with our results that more satellites 
were quenched as central galaxies prior to infall at higher 
mass. Adding a large sample of stellar ages (e.g., Gallazzi 
et al. 2005) directly to our group catalog would provide ad- 
ditional constraints on our derived star formation histories, 
as we will examine in future work. 
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APPENDIX A: IMPLICATIONS OF SATELLITE 
STELLAR MASS GROWTH FOR SUBHALO 
ABUNDANCE MATCHING 

In §5.2, we showed that stellar mass growth via star for- 
mation is nearly identical in satellite and central galaxies. 
We now explore the implications of this result for subhalo 
abundance matching (SHAM), and we discuss other pro- 
cesses that influence the evolution of satellite stellar mass. 

As described in §2.4, the SHAM technique for assign- 
ing stellar mass to subhalos has been successful in matching 
many observed galaxy statistics. At some level, this suc- 
cess is surprising: despite the differing physics of satellite 
and central galaxy evolution, with SHAM one assigns stel- 
lar mass to both a central and satellite subhalo based simply 
on the maximum subhalo mass (or circular velocity) that it 
experienced, M max , regardless of whether it is a central or 
satellite subhalo or its time since infall. Our result, that 
satellite and central galaxies grow in stellar mass by essen- 
tially the same amount, on average, justifies assigning stellar 
mass to both satellite and central subhalos under a single, 
simple prescription. 
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Figure Al. Fraction by which a central subhalo's -M max (top) 
and M star (bottom) at z = is higher than that of a satellite, 
for subhalos that had the same mass prior to the satellite's infall, 
according to SHAM. Solid curve shows median value and shaded 
region shows 68% distribution interval. A satellite typically has 
10 - 20% lower M max at z = than it would have if it had 
remained a central, with considerable scatter. At high mass, M a t ar 
growth is largely unaffected, but at low mass a satellite's M atar 
is systematically 40% lower than for its central counterparts. 



However, our results suggest tension with the way that 
SHAM typically is implemented, through the use of M max 
(or maximum circular velocity). Because satellite subhalos 
are stripped of their mass after falling into a host halo, their 
A/m ax almost always occurs prior to infall and remains con- 
stant thereafter. (In our tracking scheme, a satellite can grow 
in M max if it merges with another satellite.) By contrast, a 
central subhalo's M max continues to grow as its halo grows. 
Thus, if one uses the same instantaneous Af sta r — M max re- 
lation to assign stellar mass to both satellite and central 
subhalos, which is how SHAM typically is applied, this nec- 
essarily implies that satellites have grown less in stellar mass 
than central subhalos. Instead, the results of §5.2 suggest 
that satellite subhalos should have a higher instantaneous 
M 3 tar for their given M max than central subhalos. 

Under the assumption (for now) that satellite and cen- 
tral galaxies grow in stellar mass by the same amount, on 
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average, we quantify this possible tension by examining the 
differential M max growth of satellite versus central subha- 
los. For each satellite at z = 0, we record the time that 
it reached M max as a central subhalo, which typically was 
prior to first infall. We then identify all central subhalos at 
z = (discounting ejected satellites) that experienced the 
same M max at the same time, and we compute the median 
difference in M max at z — between the satellite and these 
central subhalos, which we call AM max . We can then treat 
satellite and central subhalo stellar mass growths as self- 
consistent and identical by increasing each satellite subhalo's 
M max at z = by its AM max . Using these increased M max 
values for satellites, we reassign stellar mass via SHAM to 
the whole subhalo population. The difference between these 
stellar masses and the stellar masses from standard M max , 
which we call AM star , indicates the amount by which the 
standard SHAM implementation underestimates satellites' 
stellar mass. 

Fig. Af shows AM max /M max and AM star /M sta r as a 
function of satellites' mass at z = 0. (For this exercise, we 
extend to masses below our SDSS group catalog limit). The 
top panel shows the fraction by which a central subhalo's 
M max at z = is higher than that of a satellite, for subhalos 
that had the same M max prior to the satellite's infall. At 
low mass, a satellite's M max at z — is typically 20% lower 
than its central subhalo counterparts, though with scatter 
to much higher values. This fraction declines with increasing 
M max , but only weakly because of the competing effects that 
more massive central subhalos grow in mass more rapidly, 
but more massive satellites fell in more recently. 

Fig. Af (bottom) shows the same, but for stellar mass, 
across a range that corresponds to the top panel. The 
mass dependence is now stronger, because of the chang- 
ing steepness in the M star — M max relation with mass: at 
M max > f0 12 Mq, M s tar increases slowly with M max , so a 
given increase in M max leads to a much smaller increase in 
M star , while at lower mass the opposite trend occurs. The 
underestimation of satellite stellar mass is negligible at high 
mass, but at M sta r < 1 9 M it is as high as 40%, with scat- 
ter to much higher values. Thus, to the extent that satellite 
and central galaxies mass growths are in fact identical, the 
standard implementation of SHAM significantly underesti- 
mates the stellar mass of low-mass satellites. 

Pasquali et al. (2010) also examined the differential stel- 
lar mass growths of central and satellite galaxies since infall, 
though in the context of the semi-analytic model of Wang 
ct al. (2008). They also saw a negligible difference at high 
mass and an increasing difference at lower mass. However, 
their difference was driven primarily by the rapid quench- 
ing of satellites after infall in their model, which we have 
argued cannot match the observed SSFR distribution (see 
also Weinmann et al. 2006; Font et al. 2008). Our empirically 
motivated approach indicates that the stellar mass growths 
via star formation should be nearly identical, on average, 
implying that satellite and central subhalos lie on separate 
instantaneous M s tar — M max relations, a scenario that can, in 
principle, still match spatial clustering measurements (Neis- 
tein et al. 2011). Promisingly, Yang et al. (2012) recently 
implemented a self-consistent approach to SHAM, in which 
they found evidence that satellite and central subhalo stellar 
mass growths are similar, in support of our results. 

However, an important caveat to Fig. Af is that star 



formation is not the only process that affects satellite stel- 
lar mass evolution. Satellites can merge with one another 
(Angulo et al. 2009; Kim et al. 2009; Wetzel et al. 2009b,a), 
though we expect that the impact of satellite mergers on 
stellar mass growth is subdominant at these masses (Wetzel 
et al., in prep.). More importantly, satellites can lose stellar 
mass via tidal stripping. As outlined in §2.3, we use a sim- 
ple binary procedure for merging/ disrupting highly stripped 
satellite subhalos, but we do not account for partial strip- 
ping of stellar mass. There are many reasons to expect that 
surviving satellites have been at least partially stripped, 
based on both the properties of satellites (e.g., Yang et al. 
2009; Kim et al. 2009; Pasquali et al. 20f0; Simha et al. 
20f0) and the presence of diffuse intracluster light (ICL) 
from stripped satellites (e.g., Willman et al. 2004; Conroy 
et al. 2007; Rudick et al. 2009; Puchwein et al. 2010). ft 
is possible that the canceling effects of stellar mass growth 
via star formation and stellar mass loss via tidal stripping 
cause satellite and central subhalos to lie on essentially the 
same instantaneous M st ar — M max relation (see also Simha 
et al. 20f0). If true, then SHAM ascribes the correct stel- 
lar masses to satellites though a fortuitous coincidence: the 
retarded growth of a satellite's M max , and thus its M sta r, 
accurately captures the physical process of tidal stripping. 
If so, then Fig. Al must also indicate the average amount of 
stellar mass stripping that surviving satellites at z = have 
suffered. Direct observations of the stellar mass fraction in 
ICL are ~ 10-25% (Zibetti et al. 2005; Gonzalez et al. 2007; 
Krick & Bernstein 2007), broadly consistent with Fig. Al. 

Finally, we note the implications of these results to 
SHAM as applied to dwarf satellites in the Local Group. 
While Fig. Af extends only down to about the mass of 
the Small Magellanic Cloud, it also indicates the increasing 
importance of satellite stellar mass growth at lower mass. 
Thus, understanding the relative importance of stellar mass 
growth and tidal stripping becomes even more important 
for dwarf satellites in the Local Group. If local dwarf satel- 
lites have not been stripped of stellar mass significantly, as 
might be implied by their tight luminosity-metallicity rela- 
tion (Kirby et al. 2011), then Fig. Al implies that the usual 
SHAM approach may underestimate their stellar masses. 
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